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Abstract 

In constructing collinear operators, which describe the production of energetic jets or energetic 
hadrons, important constraints are provided by reparametrization invariance (RPI). RPI encodes 
Lorentz invariance in a power expansion about a collinear direction, and connects the Wilson 
coefficients of operators at different orders in this expansion to all orders in as- We construct 
reparametrization invariant collinear objects. The expansion of operators built from these ob- 
jects provides an efficient way of deriving RPI relations and finding a minimal basis of operators, 
particularly when one has an observable with multiple collinear directions and/or soft particles. 
Complete basis of operators are constructed for pure glue currents at twist-4, and for operators 
with multiple collinear directions, including those appearing in e~^e~ — > 3jets, and for pp — > 2jets 
initiated via gluon-fusion. 
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I. INTRODUCTION 



Factorization theorems play a crucial role in our understanding of QCD [l|, llj. For pro- 
cesses with large momentum transfer or energy release they provide a separation of the 
high-energy perturbative contributions from the low energy process independent functions 
describing non-perturbative dynamics. The soft-collinear effective theory (SCET) provides 
a systematic approach to the separation of hard, soft, and collinear dynamics in processes 
with energetic hadrons or jets 

aaaa. 

It has an operator based approach to hard-collinear 
factorization which provides a simple framework for deriving the convolution formulae con- 
necting Wilson coefficients and collinear operators. The hard Wilson coefficients describe 
the short distance process dependent contributions, and the operators built out of collinear 
and soft fields encode the longer distance hadronization into individual energetic hadrons, 
energetic jets, or hadrons with soft momenta. With more than one collinear direction the 
factorization for SCET operators was first considered in Ref. 0], and it was demonstrated 
that the leading order operators efficiently encode traditional factorization theorems for pro- 
cesses like Deep-Inelastic Scattering (DIS), Drell-Yan, Deeply- Virtual Compton Scattering 
(DVCS), and exclusive form factors with hard momentum transfer. Compared to more tra- 
ditional methods, an advantage of the effective theory approach to high-energy factorization 
is the systematic description of power corrections by higher order operators and effective 
Lagrangians jl, [o], [l^, 11 . 



An important constraint on the construction of both leading and power suppressed 
operators in SCET is provided by reparametrization invariance (RPI). The utility of 
reparametrization invariance was first discussed in Ref. 12|] in the context of heavy quark 
effective theory (HQET). In HQET there are 3 generators for RPI, and the transformations 
involve a time-like vector where f ^ = 1. For collinear operators in SCET, RPI transfor- 
mations act on null vectors and where n-n = 2 and there are 5 generators for each 
type of collinear field. Reparametrization invariance in SCET was first discussed in Ref. 



and generalized to the complete set of RPI transformations in Ref. 13 



To see how reparametrization constraints come about, let's consider a process with mul- 
tiple energetic jets defined by an infrared safe jet algorithm, as pictured in Fig. [H We 
assign labels n^, ni^, rig, ... to the jets, which are null nf = 0, and whose vector compo- 
nents identify the directions fii of the total momentum vector of all hadrons in the jet. The 
hadronization in each jet takes place in a collinear cone about each nj, and we refer to the 
energetic particles in this jet as nj-collinear. Interactions between particles in different jets 
can take place only by hard exchange at short distance or by soft exchange at long distance. 
The description of the physics of a jet is simplified by a suitable set of coordinates, which 
are provided by nf, and a complementary null vector nf where n"^ = and rii-Ui = 2. The 
momentum of a particle in the i'th jet can be decomposed in these coordinates as 

P'' = n^■p^+nrp^+p^^^^. (1) 

The collinear modes for the jet have momentum scaling as (n^ ■ p,ni ■ p,Pni±) ~ <5(A^, 1, A) 
where A -C 1 and is a large perturbative momentum scale (the jet-energy). In cases 
where our discussion is generic to any one jet we will leave off the subscript i, so — »• n 
and Hi — > n. The definition of _L in Eq. ([1]) is relative to n and n, and for this reason we use 
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FIG. 1: Three collinear jets labeled by vectors n^. 



the notation with n-pni_ = n-Pn± = 0. When it is clear which n and n we are referring 
to we will sometimes write for p'^j_. For each n, collinear operators are built up from 
quark ^„ and gluon fields, which are labeled by their collinear direction, and describe 
quantum fluctuations close to the direction n with offshellness <^ Q^. Two collinear 
directions are described by distinct collinear fields when rij • rij ^ for i ^ j |[7|]. In Eq. ([1]) 
fii is introduced solely to provide a basis vector for the decomposition, unlike rij which has 
a physical association. For multiple collinear directions we have the freedom to introduce 
multiple raf vectors. 

Reparametrization constraints arise because the decomposition in Eq. ([T]) is not unique. 
We can shift Ui by a small amount and still have a suitable basis vector for the i'th jet. 
We also have a large amount of freedom in the choice of n^. For each {n,n} pair the most 
general set of RPI transformations which preserves the relations = 0, = 0, and n-n = 2 
are 

where the five infinitesimal parameters are {A^, e^, a}, and satisfy n-e^ = n-e-^ = n-A-^ = 
n-A"*- = 0. To ensure that n provides an equivalent physical description of the collinear 
direction for these particles requires the power counting {A^,e^,a} ~ {A^,A°,A°} (isf . 
Thus n can only be shifted by a small amount, while parametrically large values of a 
and ej[ are allowed. In B-meson decays, constraints from reparametrization invariance in 
SCET have been derived for heavy-to-light currents with parameters v and n, at the first 
subleading order in Refs. j^, 0, E3], and to second order in Ref. 15]. Results for light- 



light SCET currents with one collinear direction n, were derived at first subleading order in 



Ref. 161]. The extension of RPI relations to collinear operators involving light quark masses 



was developed in Ref. [17 



The goal of our paper is to provide a simple procedure for constructing the RPI- 
completion of operators 0{ni,ni,Vi) that depend on multiple light-like vectors {nj,nj} and 
time-like vectors fj. The procedure should be sufficiently general to be used for any hard- 
scattering process, and also easy to extend to any desired order in the twist or A expansion. 
To achieve this we must deal with a technical obstacle: so far all applications of RPI to 
hard-scattering in the SCET and in other factorization literature have constructed a com- 
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plete basis of operators first and then dealt with deriving connections between the operators 
order by order in the A expansion. This approach quickly becomes cumbersome at higher 
orders or when dealing with operators with multiple directions. For example, in this ap- 
proach the RPI completion of a basis of three jet operators 0{ni,n2,ns,ni,n2,n3), would 
require studying three copies of Eq. ([2]) or nine transformations.^ 

For cases with multiple time-like vectors, an alternative approach is known from 



HQET 131]. Here an RPI heavy quark field is constructed at the beginning, H^, which 
has an expansion that starts with the standard HQET field, Hy = + . . .. A basis of 
reparametrization invariant operators built from automatically encodes the RPI rela- 
tions at any order in the power expansion, and when expanded generates a series of opera- 
tors with connected Wilson coefficients. In this paper we develop a suitable set of RPI and 
gauge invariant objects for SCET. These objects include a quark field operator \E'„, a gluon 
field strength operator Q^'^, and 5- function operators which pick out the large momenta of 
collinear fields. The gauge invariance of these objects is ensured using a "reparametrization 
invariant Wilson line" operator VV„. These objects allow us to extend the invariant operator 
procedure to processes that depend on null-vectors. 

In hard-scattering processes, DIS provides a familiar context where the construction of 
a minimal operator basis requires judicial use of the quark and gluon equations of motion. 



and an invariance under reparametrizations of a light-like direction 19|, |20 



review see 



24 



21 



22 



23 , for a 



The invariance under reparametrizations becomes more valuable at higher 
orders in the expansion, being particularly constraining on the basis of twist-4 operators 
derived in 
Refs. 



20l . |21| . |22| . |23| |. We derive RPI constraints for collinear operators in DIS and 



compare to these classic results as a test of our setup. For DIS the minimization of the basis 



of RPI operators is quite similar to the reduction of operators in Ref. 2]J. On the other 



hand the basis of SCET operators are comprised entirely of analogs of "good" quark and 
gluon fields, namely a two-component quark field Xn and just two components of the gluon 
field strength, B!^j_. These objects both incorporate Wilson lines, and for these operators it 
is easier to find a minimal basis. The RPI relations provide Lorentz invariance connections 
between the Wilson coefficients in this basis. These constraints carry a process independence, 
they depend on the type of operators being considered, but not on the precise process in 
which they will be used. It should be emphasized that when matrix elements are considered 
for a particular process, a further reduction in the number of independent hadronic functions 
becomes possible. For twist-4 quark operators in DIS this type of further reduction was 



25|, but this type of reduction 



discussed in detail in Ref. 23|] and for inclusive B-decays in. 
is not our focus here. 

Our construction is general enough that it applies not just to DIS like processes, but 



In Ref. 18| it was shown that the construction of heavy-to-light operators can be simphfied if only 
operators in a particular frame are required, by taking linear combinations of the RPI transformations 
that only act in this frame. In Ref. |15| this was described as the derivation of RPI conditions on a 
projected surface, and the complete set of such transformations was used for the 0{X^) analysis done 
there. The formalism derived here makes a full analysis sufficiently simple that the consideration of 
projected surfaces becomes unnecessary. 
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to operators with multiple collinear directions, which are useful for processes with multiple 
hadrons and jets. These operator bases provide a starting point for deriving appropriate 
factorization theorems for different processes. The invariant operator procedure becomes 
more and more efficient as the number of directions grows. 

The outline of our paper is as follows. In section [TTl we review ingredients from SCET 
needed for our analysis. We divide hard interactions into two categories, those with an 
external hard leptonic reference vector q^, and those where the hard interaction is between 
strongly interacting particles. Since most SCET applications focus on the former case, we 
address some of the additional notational complications that occur for the latter. Section 
III introduces a formalism for using reparametrization invariant objects in SCET. A set of 
RPI invariant collinear objects is constructed in section IIII Al followed by a summary of 
identities that can be used to reduce the operator basis in section IIIIBI The inclusion of 
mass effects is considered in section IIII C\ and the expansion of the RPI objects is carried 
out in section UlI D[ Applications for constructing operators are considered in section HVl In 
section HVAI we verify that our approach provides a simple way to reproduce the known RPI 
result for the chiral-even scalar current given in Ref. 16j. In section HVBI we construct a 
general basis of field structures involving up to four active quark or gluon operators, and with 
up to four distinct collinear directions. In section HVCI we consider the special case of quark 
operators for DIS at twist-4 with one collinear direction, and compare with the literature. 
In section HVDI we derive a basis of operators for pure gluon scattering in DIS up to twist-4. 
Finally we apply the formalism to jet production. In section HVEI we demonstrate that very 
little information is gained about the operator basis describing e^e~ 2 jets. In section lTVFI 
we show that RPI turns out to be quite powerful for constraining the e+e^ 3 jet operators. 
Finally we show that RPI is also useful for two jet production from gluon-fusion, gg qq, 
and we construct a basis of operators for this process in section IIV Gl Conclusions are given 
in section IVl 



II. REVIEW OF SCET 

In sections III Al and IIIBI below, we introduce some basic definitions and properties of 
SCET that we will need for our computations. In section III CI a brief review of the null 
and time-like RPI transformations is given, and in section III Dl a review of hard-collinear 
convolutions is given since they play an important role in subsequent sections. 



A. Fields, Wilson lines, and Power Counting 

SCET fields include collinear gluons Ai^ and collinear quarks ^„ for each distinct direction 
n. An important attribute of the collinear fields is that they carry both a large label 
momentum p and a coordinate x, such as C,n,p{x)- The label momenta are picked out by 
momentum operators, Vn^n,p = n-p^n,p and 'P'^iEn,p = P±C,n.p, while derivatives id'^ act on 
the coordinate x and scale as id'^ ~ (see Ref. p]). Having two types of derivatives makes 
it simple to couple collinear and ultrasoft particles for SCETj, including ultrasoft gluons At^^ 
and quarks and when appropriate, heavy quarks /i""* as well. The soft fields for SCETn 
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are A'^, g^, and a heavy quark h^. 

We define coUinear covariant derivatives as 

in-D^ = Vn + gn-An,p, iD^^ = V^^ + g A^i^ , in- = in-O + gn- An,p . (3) 

When integrating out hard offshell fluctuations and constructing gauge invariant structures 
in SCET, it is necessary to include colhnear Wilson hues, Wn, deflned by 

Wn{x) = [ 5Z exp n-An,p{x) )] . (4) 

perms 



The colhnear flelds are deflned with the zero-bin procedure 26|. To couple ultrasoft 
degrees of freedom we deflne an ultrasoft covariant derivative 

= td^^ + gA>:^, , (5) 

that can act on colhnear flelds. At lowest order the coupling to n-colhnear flelds involves 
n-Dus and can be removed from the Lagrangian by the BPS fleld redeflnition [gI 

^n,p{x) — > Yn{x)^n,pix) , A^^g^x) Yn{x) An^q{x)Y^ (x) , (6) 

with the ultrasoft Wilson line 

Yn{x^) = Pexp(ig j ds n-Ausix^" + sn^")) . (7) 

This fleld redeflnition allows us to organize power corrections as gauge invariant products of 
colhnear and ultrasoft flelds as we discuss in the next section. In describing e^e~ — > jets in 
SCET it is convenient to make a fleld redeflnition with Wilson lines over (0, oo) rather than 
the Yn shown in Eq. ([7]) [27|, l28|. For qq jets one can use lines over (0, oo) or (— oo, 0). 
The flnal results are always independent of the choice of the reference point for Y in the 
fleld redeflnition (the — oo in Eq. (JTj)) since it does not dictate the direction of the lines in 
the flnal result (though the same choice should be used in all parts of the computation). 

Operators are formed from products of the above flelds, and the power counting for an 
operator is determined by adding up contributions from its constituents. The power counting 
for the flelds and derivatives in SCETj is^ 

^n-A, (n-A„,n-A„,A^) ~ (A2,1,A), ~ Z^"' ~ ^^ ~ A^ , 

^a'^~A^ (2n-a,n-P,P,^) ~ (A^1,A), ~ y„ ~ A° . (8) 

Here the ultrasoft flelds describe fluctuations with offshellness much less than the colhnear 
particles. These objects can be used to construct operators for processes with multiple jets. 
For a colhnear jet we have A ~ A/Q with Aqcd <C A <C Q. For a colhnear hadron we have 
a smaller A, namely A ~ Aqcd/<5- For processes with two or more hadrons the interactions 



We will often suppress the labels on coUinear fields when writing them out is not essential. 
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in the theory SCETn must be considered. With a small parameter 77 ~ Aqcd/Q ^ 1 the 
power counting of fields in this theory are 

«9,^~77, (27i-9,7i-P,P„^)~ (772,1,77), W^r^S^r^rf. (9) 

Here the soft fields describe fluctuations with similar offshellness to the coUinear fields. In 
cases with jets and energetic hadrons a succession of SCETi and SCETn theories needs to 
be considered. 

Our article focuses on building reparametrization invariant operators from products of 
collinear fields that describe an underlying hard interaction, since this is the most involved 
part of the construction. The simple strategy we follow to incorporate "ultrasoft" and "soft" 
fields into the analysis is summarized in sections III Bl and III CI below. 



B. Gauge Invariant Field Products and Convolutions 



To build operators in SCET we want to use structures which are gauge invariant and 
homogeneous in the power counting. Although the precise manner in which the Wilson 
lines Wn appears is determined by matching, and the precise manner in which Wilson lines 
Yn appear is determined by ultrasoft-coUinear factorization, some general structures can be 
identified. For SCETj a convenient set of structures are: 



(10) 



together with the label momentum operator and derivative operator id^ acting on these 
gauge invariant structures. The collinear fields in Eq. ffTOl) are the ones obtained after the 
field redefinition in Eq. It is convenient to be able to switch the collinear derivatives 
multiplied by Wilson lines for gauge invariant field strengths, for which we use 



in-'Dn = in-d + gn-Bn 



in-T>n = in- d — gn-Bn , 



and note that n-Vn = Vn- Here the field strength tensors are 



In 



gn-Bn 



=-[in-Vn,in-Vr, 



(11) 



(12) 



where the label operators and derivatives act only on fields inside the outer square brackets, 
and gB'^j_ and gn-Bn are Hermit ian. 

For SCETn with hadrons we have the same collinear invariant objects as in Eq. fl2^ . and 
similar soft invariant objects, that are obtained by replacing the ultrasoft fields by their soft 
counterparts, K = (F>r) - K = (SlK), V^^ P^, and q^, = (F^t^,^) ^ ^ ^^t^^). 
The soft Wilson fine S"^ is generated by integrating out offshell fluctuations which determine 
its direction n, and outgoing/incoming boundary conditions. Most often these operators 
can be constructed by a matching calculation from SCETj, in which case the properties 



of the soft Wilson lines are directly inherited from the ultrasoft ones in SCETi 29[, and 
the product of C{Q^,uJi) from Eq. (|20|) and J{ui,kj) from Eq. fl28l) becomes the Wilson 
coefficient of the factorized operator in SCETn. In this paper we focus on SCETi examples. 
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n 
n 



n-Dr, 



n-Dr, 



in 

w 



Type (I) 



n + A-L 
n 



nl. 



— n-Dn - 



n-Dn 

w 



■D„ 



n 
n 

n-Dn 

D^. 

n± 

n-Dn 

in 

w 



Type (II) 



n 

n + £-* 
n-Dn 



n-Dn'—e^-D. 
n-Dn + e^-D^ 

1 + ^^^^^^)^" 



n-Dr, 



--D^)W 



TABLE I: Summary of infinitesimal type I and II transformations from Ref. 13] • With multiple 
collinear directions these transformations exist for each {nj,nj} pair. 



C. Reparametrization invariance 



When a set of fields have their largest momentum component in a light-like or time-like di- 
rection then the structure of operators built from these fields is constrained by reparametriza- 
tion invariance. This invariance appears due to the ambiguity in the decomposition of 
momenta in terms of basis vectors and in terms of large and small components. For a 
collinear momentum, the set of five transformations on the light-like basis vectors raf and 

were given in Eq. ([2]). These infinitesimal changes preserve the relations nf = 0, nf = 0, 
Hi- Hi = 2, and with the power counting {A-^,e^,a} ~ {A, A°, A°} can have no physical 
consequences on the description of an observable. The type-Ill boost simply ensures that 
{jj^Nui) — {jj^Nfii) — {jj^Drii) + {^Dfii) = for each i, where {jj^Nn.i) counts the number of 
Hi factors in the numerator of an operator, {^Dfii) counts the numbers of fij factors in the 
denominator, etc. With three collinear directions an example of a type-Ill RPI invariant 
parameter is 

The type-I and type-II transformations of collinear objects are more interesting and are 
summarized in Table [H which we take from Ref. . Since the factors induced by these 
transformations occur at different orders in A, demanding overall invariance of a physical 
process provides connections between the Wilson coefficients of operators at different orders 
in the expansion. 

When we couple collinear and ultrasoft particles there is another ambiguity, associated 
with the decomposition of a collinear momentum into large and small pieces. If the total 
momentum of a collinear particle is decomposed into the sum of a large collinear and 
a small ultrasoft momentum k^: 



n 



—n-{p + k) + -i^n-k + {p± + fcj 



(14) 



then operators must be invariant under a transformation that takes n-p n-p + n-i, 

n-k — n-i, and k'^ k^ — To construct invariant objects that have 



Pi 



p1 + ^ 



n-k 
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nice gauge transformation properties we use the combined covariant derivatives [1 ll. |30[| . 

^D^^ + WniD'^^^^Wi , in-D^ + Wntn-DusWl . (15) 

This can be implemented by taking 

zP;:^ ^ tDf^, = tV^^ + tDi^ , V„. ^ m-DMi = P„ + tn-D^, , (16) 

and then expanding in A. The results in Eq. (fT6l) give powerful relations as they relate the 
coefficients of operators involving coUinear fields to those involving ultrasoft fields. These 
relations are quite easy to derive order by order in A. Note that reparametrization constraints 
associated with transformation of the ultrasoft Wilson line ¥„, are automatically enforced by 
the other constraints.^ 



Finally we review RPI for a time-like vector from HQET [1^. The momentum of a 
heavy quark is decomposed as = mv'^ + k'^, where m is the heavy quark's mass, v'^ is its 
velocity, and k^^ is a residual momentum of order mA^. For an infinitesimal (3^ ~ A^ with 
v-P = 0, the shifts 

v^^v^ + f]^ and k^^k^-mf^", (17) 

can have no physical consequences. This implies invariance under the infinitesimal change 
hy ^ hy + 6hy with 6h^ = {im(3-x + $j2)h^. A superfield can be constructed which is 
invariant under the full transformation Il2| 



H^{x) = e" 

where 



■imv-x 



1 ll\ 



v^(l + i;-V/|V|) V \V\J 



1 + 



Kix) (18) 



= t;^ + iD'^Jm . (19) 



Using this superfield one can build operators O = 0[Hv{x),Df^] that are invariant under 
reparametrizations of the time-like vector. Here = e~*"*^'^[l + ip/[2m) + . . .]hy at the 
first non-trivial order. Note that for heavy quarks, no dynamic component of the momentum 
is the same size as the hard fiuctuations, so there is no analog of the 5-functions in Eq. fl2^ . 
This is the main complication we face in constructing invariant operators in SCET. The 
closest one gets in HQET is when we have two auxiliary time-like vectors, v and v', such as 
in P ^ D^*^ decays. Here the invariant Wilson coefficients must be functions C(V-V') jsij. 



D. Convolutions 

In the presence of collinear fields a hard interaction can introduce convolutions in variables 
Ui between the perturbatively calculable Wilson coefficient C{Q'^, uoi) and the matrix element 



^ For example, prior to the field redefinition only the combination in-T) = in-d + gn-Aus + gn-An appears 
acting on collinear fields. A type-I transformation connects this to a Pj;, and Eq. (|16p then connects this 
to the same iD:^^ that one would find by direct transformation of K„. 
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of the collinear operators. In this case the amphtude, cross-section, or decay rate has the 
form 



A 



[duji ■ ■ ■ duk] C{Q^, uji) {0{uJi)) . 



(20) 



The convolutions occur because a component of the hard momentum and of one or more 
colhnear momenta are (9(A°). The exchange of momentum between the hard and colhnear 
components yields a convolution in variables Ui, where the number of such variables is 
constrained by gauge invariance and by momentum conservation in the matrix element. A 
gauge invariant momentum from the collinear fields can be picked out by a delta function 
acting on one of the collinear objects in Eq. ffTOl) . such as [6{u — n-Vn)Xn], and traditionally 
in SCET a subscript notation is used for these products, 



Xn,L 



S{uj- Vn)Xn 



{gn-Bn)uj= gn-Bn6{uj -V\) 



(21) 



We will refer to these as homogeneous objects since they have a definite order in A, and call 
the operators build from these objects homogeneous operators. As an example we have the 
bilinear scalar operator. 



0{Ui,UJ2) 



Xn,u)i 2 Xn,uJ2 ■ 



(22) 



When we consider RPI it will be convenient to use different 5 functions and convolution 
variables u, that are type-III invariant. Essentially each Vn = n-Vn must be multiplied by 
a scalar transforming as n under RPI type-III. There are two cases to consider: 



i) situations where there is a reference vector q'^ for the hard interaction, 
Aqcd; which is external to the QCD dynamics. 



» 



ii) situations where the hard interactions are purely from strongly interacting particles. 

Case i) applies to examples such as DIS where q'^ is the momentum transfer from the 
virtual photon, or e~^e~ — >■ jets where q'^ is the four momentum of the e~^e~ pair. Here 
we can use n-q ~ A*^ to make the 5-function type-III invariant for n-collinear fields. Since 

^ AAqcd ^ ^QCD ^6 know that n-q ^ n-p, where p is the momentum of a collinear 
particle in the jet. Thus we use a variable a; with mass dimension two, and will find 6- 
functions of the form^ 



6[uj - n-qVn) 



(23) 



^ For _B-decays these type-III invariant (5-functions were used in Ref. [14|, with ~ rribv'^'-, 6{uj — n-qVn) 
S{uj — mbfi-vVn) — l/mbS{Lu' — n-vVn), where Cj — miyUj' . This form of invariant (5-function was also 



quite useful for analyzing the factorization theorem for 



J/^jX in Ref. 32]. 
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We also introduce a subscript notation with hatted variables, 



Xr. 



i9B: 



5{LJ-n-qVn)Xn 



(gn-Bn) 



gn-BnS(uj-n-qVl) 



(24) 



Since (5(a) — n ■ qVn) ~ A'^, it is leading order in the power counting. Furthermore, we have 
6{iJj — n ■ qV) = 5{Cj/n ■ q — V) /n ■ q, so identifying ci) = n-q uj there is no real change to 
the structure of Eq. ( !20l) . An operator built out of the components given in Eq. ( !24l) has 
multiple labels, 0{uJi,Co2-, ■ ■ ■)■, and the Wilson coefficient for the operator will be a function 
of the same parameters, C(ci)i,c2'2, • • ■)■, yielding Eq. (120|) with a;'s replacing cj's. 

For processes in case ii) there is no analog of the external q^. Examples here include 
PP — > jets, or any other hard process that does not involve external leptons or photons. 
The key difference with case i) is that here the hard interaction must involve two or more 



coUinear directions, so we are guaranteed that there are scalar products rii 



A°. For 



this type of reaction the type-III invariant (5-functions which are convoluted with Wilson 
coefficients always involve large momenta for two different collinear directions. 



(25) 



Here P„. acts on a gauge invariant block of nj-collinear fields, and Vn^ acts on a block of 
nj-collinear fields. Since this 5-operator does not act on a single block of collinear fields 
we will not use a subscript notation like Eq. (1241) for Uij. In this case the structure of the 
factorization theorem between operators and Wilson coefficients is a bit different than in 
Eq. fl20|) . For example, consider an operator with collinear objects for four directions, where 
the convolution is 



[H du,,] C{u,,)[l[ A, J Xn^gBiJigB, 



n4lXn2 



(26) 



km 



Here the products are over the six unique pairs ij with i ^ and in the A^m acts on 
the nj-coUinear field(s). The convolutions in Eq. fl26|) can be manipulated into the form of 
Eq. ( !20!) by inserting four factors of 1 = / duJidiuJi — Vm), writing 5^- = (5(a)jj 
and carrying out the integrals over the six ujij^s to give 



rii-rij uJiUJj/2) 



j [dui ■ ■ ■ dui] C[ni-nj i^ii^j)Xm,u^AgBi,^^^){gB^^ ,^^ 



)X 



(27) 



Here the RPI-III transformation of the measure cancels against that of the (5-functions in 
the operator, and RPI has constrained the Wilson coefficients to only depend on invariant 
products ?T,i-?T,2Ci;iCi;2, ni-n^uJiLJs, etc. 

Due to the simplicity of the ultrasoft-collinear coupling at leading order in SCET a further 
factorization of the EFT matrix element can be made into collinear pieces J, and ultrasoft 
pieces S at each order in the power counting: 



(0(0.,)) 



(28) 



However it is the factorization in Eq. fl20l) that will be central to our discussion of 
reparametrization invariant operators. 
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III. REPARAMETRIZATION INVARIANT OBJECTS FOR SCET 



To construct an expansion in operators in SCET the standard procedure is to build a 
gauge invariant basis of operators with definite power counting, order A*', and to assign a 
Wilson coefficient to each one. Afterwards one can impose RPI order by order and find 
relations among Wilson coefficients. On the contrary what we will do is to start with RPI 
and gauge invariant objects, to be constructed in section [Til A[ These objects do not have a 
definite power counting order, in particular we will know the order in the A-expansion where 
they start, but they will contain terms at all higher orders as well. We build a basis with 
these RPI and gauge invariant objects, which is made minimal using equations of motion and 
kinematic constraints as discussed below in section IIII B[ (Equation of motion constraints 
for homogeneous operators are also summarized in this section.) Each element of this basis 
is assigned a Wilson coefficient, and then the elements are expanded to find the final basis 
with elements of a definite power counting. In this way we immediately obtain relations 
between Wilson coefficients of operators at different orders. Once we expand and check 
for redundancy, the number of independent Wilson coefficients is equal to the number of 
independent RPI operators in the reduced basis. 



A. Construction of RPI and Gauge Invariant objects 



We now construct reparametrization invariant objects in SCET whose leading terms give 
the fields in Eq. (jH]). These are then generalized to objects that are simultaneously RPI 
and gauge invariant whose leading terms give the objects in Eqs. f fTOfMl) . For simplicity 
only coUinear objects are considered in this section. Pulling out the large phases from the 
collinear quark field and gluon field strength, and decomposing the full theory field into 
independent collinear sectors we have at tree level, 

^(x) = ^ e-*"-^"7/'„(x) , C'^ix) = J2 e"*"'^"G'r(a;) . (29) 

n n 

Full Lorentz invariance act on the fields il){x) and G^'^{x), but the RPI transformations that 
we are interested acts independently on each collinear sector labeled by n. Two sectors i, j 
are independent if rii ■ rij ^ A^, and the sums in Eq. fl29|) are really over equivalence classes, 
{n}, where a class consists of vectors related by RPI. From the discussion in section Hi CI the 
n-reparametrization invariant collinear quark and field strength are easy to identify 

= (l + -^Pn^) , ^9 G^, ^ [iD^, iD^] . (30) 



Under the transformations in Table [T] for {n, ra}, the quark field ipn remains invariant 13 
while the gluon tensor is invariant because the vector is invariant. To make the fields in 
Eq. (l30l) invariant under the additional reparametrization transformations that link collinear 
and ultrasoft derivatives we replace in-Dn in- Dn + gn- A^s, iD'^± ^-Dn± + ^«^-^us±^n' 
and in-Dn in-Dn+Wnin-DusW\. After this replacement the decoupling field redefinitions 
in Eq. ([6]) can be made. In Eq. fl30|) n/t^n = 0, and the term in ipn with a _L-covariant derivative 
corresponds to the two components of the full fermion field that are small when p±/n-p <^ 1. 
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Since t^ipn 7^ 0, the ipn field does not provide a definite power counting for operators. For 
example, Tpn^iipn ~ A° whereas ipni^i^n ~ A"^. 

We also need reparametrization invariant 5-functions whose expansions reproduce 
Eqs. fl23l) and fl25|) at lowest order. For example, these are needed to construct an RPI 
operator which when expanded gives Xn,uj-i_^Xn,uj2 at lowest order. For situations where there 
is an external hard vector the invariant 5-function is 

A, = 5{iu - 2q ■ idnj = 6{lu - Ui-qVnJ + ■ ■ ■ , (31) 

where as described in section IIIB[ is a parameter specific to the kinematics of the process 
being studied. Notice that 6{lj — 2q ■ idn) starts at (9(A''), is RPI, and is gauge invariant 
when acting on singlet operators. Here 

= Y^n + P:^ + ym-9., (32) 

and functions of idf^ ~ (A^, 1, A) can be expanded in powers of A. Note that Vn and V^^ 
are only non-zero when they act on ri-collinear fields. It is useful to extend this property to 
the full id'^, which we can do by distributing an id^ derivative across all fields that it acts 
on, writing for example id^ijjni4'n2 = (^f^niV'ni)^n2 + i^n-i_{idn2'^n2)- some hard processes 
there is more than one external hard vector, and a natural question arises as to whether 
provides a unique choice for this construction. For example, in DVCS, 7*p we 
have the momentum q^ of the incoming 7* and the momentum q'^ of the outgoing 7*^*^ In 
Appendix |X] we show that as long as — ~ A or smaller, the choice g suffices, since for 
the purpose of constructing a basis of operators it is equivalent to the choice of any linear 
combination of g and g'. On the other hand, for situations where there is no external hard 
vector g^, the appropriate RPI 5-function is 

kij = 6{uij - 2idn, ■ idn^) = S (^Uij - ^rii ■ rij Vn{Pn^ + ... . (33) 

This 5-function operator acts on two independent collinear directions. In general we must 
include in an operator a set of Aj and Ajj which are linearly independent. Once we expand, 
the first term in the series for Ajj is not independent of the first term from Aj, so the 
(5-function shown on the RHS of Eq. fl33|) can always be eliminated, as we did in Eq. fl27|) . 

We will also make use of a reparametrization invariant Wilson line, W^, which has the 
same gauge transformation properties as Wn, 

Wn = e-*^" . (34) 

Here the operator Rn starts with a term at 0{\) and is built of n-coUinear gluon fields, 

Rn = Rn[Vn.V^n±^9KLX], (35) 

where the vector is either q^ or id'^, with n ■ n' \^ . Furthermore, i?„ is Hermitian, 
dimensionless, and collinear gauge invariant. We leave the explicit construction of Rn to 
section IIII Dl below, and for the remainder of this section take these properties as given. 

Under collinear gauge transformations, ipn and W„ transform the same way as ^„ and 
Wn, and G^*^ transforms as a nonabelian field strength. Thus using we can form analogs 
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of the results in Eq. that are simultaneously RPI and gauge invariant, namely the 
superfields 

vi/„ ^ wiiJn , = wiGrw. . (36) 

For cases with an external we also introduce a subscript notation, 

^n,^ = - 2q-td^)^n] , QZ(. = [e^r K'^ + 2q-iK)] . (37) 

Operators built out of the superfields and Q^'' are simultaneously RPI and gauge invari- 
ant. They are not homogeneous in the power counting, but the superfields reduce to the 
objects in Eq. (12^ at lowest order in the A expansion. For example, the superfield for the 
fermion 

,.H,.|.a.,8,j ^1 ^ i-,9i|) ^„ 



= Xn + • • • • (38) 



Similarly, {gyyin^)ig{Qf^''') = VnQB^y + . . .. Thus to form a RPI version of the bilinear 
fermion operator 0{ijJi,uj2) in Eq. (!22l) we simply take 

Cl{uJi,UJ2) = ^n,CoA^n,U2. (39) 

and note that expanding in A gives Q(a)i, 0)2) = {n-q)"'^ 0(^011, UJ2) + . . .. 

We will also need the equations of motion for the RPI quark and gauge superfields in 
Eq. fl5B]) . The n-collinear Lagrangian for the quark field is [i] 

C-,n = L {in ■ + ^^—i—^l^^ , (40) 

We can write Eq. fHU]) in terms of ipn as a simple Dirac Lagrangian 

Cqn = tpnipn'ipn , (41) 

The equation of motion for ipn is a simple Dirac equation pni'n = 0. Using W„W^ = 1, we 
can write yVlipnyVnyVl4'n = 0, and thus obtain the equation of motion for 

fn'^n = 0. (42) 

Here T>l^ is the RPI and gauge invariant derivative 

= W^D^ Wn = e'^" V;, e-'^" . (43) 

For the gluon field we have the equation of motion [iD^,G!^'^] = igT^J^f i'lT'^'j^ipl: 
for the superfield 

tv^g^'^ = [iv:, g:^] = -ig t^Ht^yH ■ (44) 

Note that igg^^'' = \iV^;^,iV''^]. 
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B. Reducing the Operator Basis 



In general there are three steps that one can consider to reduce the perturbative and 
nonperturbative information in the EFT to its minimal form: 

a) Find a minimal basis of homogeneous operators and of RPl operators that suffice at 
the desired order in A. The homogeneous operators can be written entirely in terms 
of x„, i3|^_L, and P^. 

b) Compare the homogeneous and RPI basis to determine which perturbative Wilson 
coefficients are fixed by RPL 

c) Consider the decomposition of matrix elements of operators in the homogeneous basis, 
and derive further relations between the resulting non-perturbative functions. 

Generically the relation between the operator basis looks like 

E E /tn d^^^ = E E /tn ^^^-i ^^(^^o ^oeiujj)] +..., (45) 

where Q^(ti)j) are RPI operators and 0£{u!j) are homogeneous operators, and the ellipse 
denotes higher order terms in the power expansion. In general our focus in this article is 
to carry out b) which is still largely process independent. For the most part we give no 
discussion of item c), which obviously must be considered process by process. In order to 
consider b) we must first determine a) which is the focus of this section. We will discuss the 
equations of motion and other relations that allow a reduction in the basis of operators at 
each order in A. 

First wc consider the gauge invariant objects with homogeneous power counting. We 
would like to demonstrate that all operators can be reduced to a form that only involves the 
basic building blocks Xn, g^n±' ^^'^ other homogeneous objects can be reduced to 

these. For example, one might think that the objects gB'^j_ = [1/V W'^[iD'^j_,iD!i!^j_]W] and 
9^±2 = [1/^ ^^[^-^n±> in-Dn\W] are independent. However they are related to the building 
blocks by 

9BT± = =Ki9B1_) - =Vl{9B'i) + p [gB'i, gBl] , (46) 
9^2 = =K{9n-B) - hn-dn{gB'J + ^[gB'^, gn-B] , 

where we will see below that n ■ B and in ■ dnB'^ can also be reduced using the gluon equation 
of motion. For Xn the equation of motion is 

in-dnXn = -{gn-Bn)Xn - ift^ift Xn , (47) 

' n 

which allows us to eliminate in-dn derivatives on Xn- To obtain the equations of motion for 
the gluon objects we consider -g^T'^J^ji^n^nT'^Wl-fy^ = [i'^uA^'^n^'^'^n]]- Expanding 
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in A and multiplying on the right with 5{uj — v\) gives three equations 

/ 

i^n■^^gB';:l = -[V^[gBl gB^]]^ - [gB^, [P^^gB^^]]^ - [gB^, [gB^,gBnl 
+ - [nM^^^^ll + ^[gB^,gn.Bl - y (gB^)^,] 



v 



1 



[m-d^gn-B]^ - [iV^fgn-B]^ - [Vi[gBl, gn-B]]^ - [gB^,, [Vlgn-B 



+ [gB^,, [in-dngBl]]^ - [gB^,, [gBlgn-B]]^ + [in-dnV^gBl]^ 



— [{gn-B)^^^', {gn-B)^] . 



(48) 



Here we sum over the color A, over the flavors /, and integrate over the repeated index 
uj'. In our analysis the flrst two equations will be used to eliminate gn-Bn and in-dngB^ 
respectively. The last relation only becomes relevant at higher orders than those we consider 
here. The above relations imply that when building a homogeneous basis of operators we 
do not need to consider the objects 



in-dnXr. 



n-B„ 



in-dnB'^j_ , 



CM 

• 



(49) 



Next we derive relations that can be used to reduce RPI operators to a minimal form. 
Given the deflnition in Eq. (143|) . we can write = id!^ + [ii'^^], and it is straightforward 
using Eq. (!66l) below to prove that 



[q-idnit>l] = Qf^igg. 



n 1 



(50) 



and hence that q^l'Dn] = 0. (The results here and below apply equally well for t = g and 
t = idn' with n ■ n' ~ A°. For simplicity we use the notation with t = q.) Eq. fl50p can be 
used to rewrite the quark superfields equation of motion in Eq. P2l) as 



Iq-idn 



QfiluigQn 



flU 



Since q ■ idn ^(t^ ~ 2g ■ idn) = 5{y> — 2g • idn) we also have the result 



q-idn'^n,c 



n,LO 



(51) 



(52) 



In a similar way, q ■ idnG: 



-Cj /'2)Qn^Cj- '^^^ coUinear gluon equation of motion for Q, 



in Eq. fl44j) can be rewritten as 



q-id, 



(53) 
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The quark and gluon operators will have u subscripts, '^n,w and ^n,tj, so only the equations 
of motion in Eqs. (1511153^ should be used to remove derivatives since the idn derivatives 
commute with the presence of the 5-function denoted by the u subscript. The QCD Bianchi 
identity, Df,G^„+D^G^f,+D^Gf,^ = 0, also gives a relation for g^", namely + 
T^n^n'^ = 0. Rearranging it gives the following relation 

(54) 

which implies that id'^Q!^^ , and id^Q'^^ are not all independent. Closing Eq. (15^ 

with 7^ allows us to remove i<^nQn'^, which is how we will choose to use this identity in quark 
operators. An analog of the Bianchi identity does not occur for the building block gBI^ in 
homogeneous operators; it easy to verify that when expanded in A, Eq. ( !54l) is trivially 
satisfied. Eqs. (l5TI - [5^ are the RPI equivalent of the results in Eqs. ( H71|48l) . and can be used 
to reduce the RPI operator basis. 

The above results imply that when building an RPI operator basis we do not need to 
consider the objects 

This list is not exhaustive. By manipulating operators in specific situations further struc- 
tures can be eliminated using a combination of the above identities. For example, for in 
sections [IV CI and II V Dl below we will see that q^Qn^id'^, with the id^ acts on a n-coUinear 
quark or gluon field, can be eliminated. 

In principle one can just count the number of RPI operators and compare to the number 
of operators in a homogeneous operator basis with definite power counting to determine 
whether there are any RPI constraints on the Wilson coefficients. The key issue here is that 
of linear independence, even if one has the the same number of operators in the RPI and 
homogeneous basis, it could be that two RPI operators constrain the same linear combination 
of operators in the homogeneous basis. 



-q-idr. 



+ 



-q-idr^ 



']q 



C. Extension to Massive Collinear Fields 



Massive collinear quarks in SCET were first studied in Refs. 33|, |3J]. After the field 
redefinition in Eq. ([6]) they have the LO Lagrangian 



qn 



in-Dn + {iIfi-m)- 



tn ■ Dr, 



(56) 



The appropriate RPI transformations with massive quarks were determined in Ref. 17 
The only change is in the type-II transformation of the fermion field, where one has to add 
a mass dependent term: 



in 



(57) 
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Under this transformation the Lagrangian in Eq. ([6]) falls into two invariant parts, one fixed 
by the leading order kinetic term and one whose coefficient encodes the choice of mass 
scheme. Note that the RPI transformation itself is not modified by the presence of a mass 
term, the transformation of n is still exactly as in Eq. ([2]). 

We can now build an analog of the RPI superfield for a massive coUinear quark. The 
reparametrization invariant quark field is 

This leads to the modified RPI superfield for a massive collinear quark 

Xn . (59) 

This result is included for completeness. Our focus in the remainder of the paper will be on 
massless collinear quark fields. 



^ r 



iRn 



I rt. ^ 



D. Determination of Rn and Expansion of '^n and Qn 

In this section we derive an expression for appearing in the RPI Wilson line, and then 
expand the invariant objects Ql^ ^ S{il! — 2q-idn), and S{uJi2 — 2i9„^ -idn.^)- We can define 
the collinear Wilson line Wn by the equation: 

[{n-D^)Wn] = Q. (60) 

We define the RPI W„ generalizing ( l60l) to a covariant derivative Dn along a (non light-like) 
direction t as: 

[(t-D„)W„] = 0, (61) 
where t is such that n ■ t ~ A". This implies the momentum space representation: 



5^ "''"((Til) 

.perms ^ ^ ' ' 



(62) 



We would like to find Rn such that Wn = WnC"*'^". Thus e~*^" is the operator that rotates 
Wn from the light-like direction n to the direction t. W„ is reparametrization invariant to the 
choice of the basis vector n, which labels the n-collinear fields A^, since such reparametriza- 
tions cannot change the fact that n ■ t ~ A°. Recall that the subscript n on >V„ labels the 
equivalence class {n} of vectors that are related by type-I and type-Ill RPI transformations. 
For any t such that n ■ t ~ A" we have 

t- An = ^n-An + ... , (63) 



t ■ idn V 

and thus 



>V„ = iy„ + . . . , (64) 
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where the elhpses represent power suppressed terms. In Eq. (163|) the n ■ t's in the numerator 
and denominator cancel out in the leading term, leaving a t independent result. 

For situations where we have an external hard vector q^, we can simply take = and 
use the corresponding Wn as the RPI invariant Wilson line. 

For situations where there is no external q'^, the choice for in Wn is less obvious since 
the only available RPI vectors are operators themselves, id'^,, where n' is a distinct collinear 
direction from n. In this situation, any choice = id'^, satisfying n ■ t = n ■ n' Vn' + . . . ~ 
A° + . . . is equally good, and the existence of the hard interaction guarantees that such 
an n' exists. In this case still yields Wn at lowest order, and hence only behaves like 
an operator in the n' direction through terms in the power corrections, namely the ellipsis 
in Eq. (l64l) . In these ellipse terms the z9„''s appear linearly order by order. Since the 
derivative idn' does not act on n-coUinear fields it behaves just like an external vector q as 
far as manipulations related to the n-collinear fields are concerned. 

In the remainder of this section we adopt the notation t = q, even though the algebra 
applies equally well to both cases mentioned above, with the substitution q ^ t = idn' in 
appropriate places. The only complication for the case t = idn' is that the dot product 
n ■ idn' must be expanded using 

fi . 77^' fi . 

2idn-idn' = -^VnVn' +n'-idniVn' +n-idn'iVn + 2idn±-idn'± H —n-idnVn' 

H — n' ■ idn'Vn + n ■ id^, n ■ idn + n' ■ idn ' ^^n' H ^— n ■ idn n' ■ idn' , (65) 

where the first term is ~ A'', the next two ~ A, the following three are ~ A^, then the next 
two are ~ A^, and the last one is ~ A^. 

Adopting t = q, Eq. f l6T]) can be used to prove that 

{q-tDn)=Wn{q-tdn)Wi. (66) 

To calculate iRn we exploit Eq. fl^B]) and calculate iRn order by order in A. Substituting 
Eq. (131 into Eq. §^ we find 

{q-tVn)=e-'''-{q-tdn)e''^. (67) 

Because of the Hermicity of iVI^ and id^, Rn is Hermitian. Applying the Hadamard formula 
to Eq. (1671) we obtain 

oo ^ 

{q ■ iVn) = {q ■ idn) + -«(^ • *^")' (^^-)' }} ' (68) 
where {{A, B}} = [A, B] and 

j 

{{A B^}} = {{[A, B], B^-'}} = [[... [A:B,]B,]...,]B] . (69) 
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Expanding i?„ in terms with Rn '' ~ A'^ we can expand all the objects in Eq. fl68|) in A and 
solve the resulting equations order by order for Rn^. Thus we write 



(g ■ iVn 
{q ■ idn 



Rn = J2^Ri'^ 

k=l 

n ■ q 



2 

n ■ q 



fl • Q Tl ' 



Vn + (g± ■ Vnl.) 



n ■ q 



n ■ idr, 



(70) 



2 ' 2 

(g ■ is a derivative operator, so when it acts in a commutator with (gB!^) we have 

[(g-^90,(^/S::)]= [g-za„(^?S::)], (71) 

where the last set of square brackets means that the derivative acts only inside. Substituting 
Eq. (!70|) into (!68|) we can solve for iR^'' . The first two terms are 



iR^''> 



+ 



n-qVj 
1 

n-qVn 
2 



(72) 



(n-g) {gn-Bn) 



4g±-Kj 
{n-qVn) 



n-qVn 



-^q^-igB^) ,q±-{gB^) 



-n-qVn 

The n ■ ;B„ term should be further reduced with the equation of motion in Eq. (l48l) to 
terms involving x„ and S^j^. In terms of the iRn'' we can determine the A expansion of the 
invariant Wilson line 



k=0 

Using the definition in Eq. the first few terms are 
The expansion of the invariant Wilson line is therefore 



Wn = Wn-WnitR^^^) + Wn 



+ ... . 



(73) 



(74) 



(75) 



Using these Rn^^^ and Table [T] it is simple to check explicitly that W„ is RPI up to order 
O(A^). Note that we did not assign a suppression for q± anywhere above (ie, we took 
gx ~ A°). Taking g^ ~ A causes further suppression of some of the terms in Eq. (!72|) . For 
cases where g_L = the expansion of W„ starts at O(A^). 

We will also need the A expansion of the invariant 5-functions, 5{uj — 2g ■ idn) and 5(0)12 — 
2idn^ ■ idnz)- For the former we have 



6{u — 2g ■ idn) = ^(^^ — n-qVn — 2g_L-'P„_L — n-qin-dn 



(76) 



n-q 



k=l 
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where the first two terms are 



p. 



(1) 



'2q±-'Pn± d 
n-q duj 



P. 



(2) = 2 



g±-Pn±\2 (p n-q d . 

^^--^nh-- (77) 

n-q J auj^ n-q auj 



When combining the operator with the Wilson coefficient C{uJi) we can integrate by parts 
to move these derivatives onto the C{u!i) and leave a simple 5-function in the operator. For 
the 5-function with two collinear directions we have 



5(a)i2 - 2idn,-idn2 



(78) 



fc=l 



1^12 



ni-n2 



2 T^n\Pn2 



where the first two terms are 



fnin2 



Vn^ni-idn2i. + 'Pn2n2-idn^_ 



d_ 

duj ' 

2 (P 



ni-n2 



J duj^ 



2idn^±-idn2± H —ni-idnj^n2 H 7^ — ^2-«<9„2P„i 



d 

- duo 



(79) 



All terms with n - idn in Eqs. ( I77|) and (17^ will be further reduced by the equations of 
motion in Eqs. fHTl) and fHHj) when they appear in operators. 
Finally we expand the superfields in Eq. fl36l) in A, writing 



(k) 



k=l 



Eel' 

k=l 



(80) 



where ^^^'1 ~ and 01^1""^ ~ The expansion of the quark superfield is straightforward, 
the first few orders are 



T{k)fj,u 



(1) 



n-q 

1 n 



(81) 



^(3)_^/-^(2) + r„(2)y 1+ \ 

~ n-q l,'^-.— + LPn Xu,.\ + ^^^^^ ^^„„ 



+ 



■q 

L 

i ^^h^ ty 

to' Z, 



1 

2 



Here there is an implicit integration over the repeated indices uja and cjf,. For the gluon 
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superfield first it is useful to expand W^Gfj,i^W: 



n, 



n„ 



[PgB^.] - ^ [PgBi] + [VgB^^^] - ^ [Pgn-B] + ^ [Pgn-B] 



4 



4 



2 

wliere gB^^ and (/^B^^ are given by tlie combinations of fields in Eq. fH6|) . Using this result 
to determine the first few terms Q^c^^ from expanding Eq. (136!) . we find 



(82) 



n,Lj 



2{n ■ q) 
1 



[n^igB: 



n' 



(83) 



n-q 



(1) 



+ ^^n^ri'^-n>^fi'^){gn-BX + | ^RZ^^.^^" (9B'nA.U-n^{gB''^^) 



^gQfr 



n'^^ujigB: 
-1 



2{n ■ q) 



n^£^uj{gB:^)^ 

[{n^Vl-n''V^^){gn-B.,)^] - [zn-d^n^gBl^ - n'' gB^^J] 

+ [n'gBl^,gn-B^] - [n'' gB^^^^, gn-B,,] - \n''u;{gB':^^)^p^^'>^ - n'^ooigB^^U^^ 
+ ... 



(2)t 



where again there is an implicit integration over Ua in terms where it appears. Here the 
ellipsis denotes terms in Qn^^'^ with an iRn'"^^ or pn^ which were not needed for our analysis. 
The (gn-Bn) and [in-dngB'^] terms are further reduced to 'Pj.'s, {gB'^Ys, and Xn's by 
using the equation of motion in Eq. ( HHi) . Finally, recall that the expansion coefficients in 
Eqs. ( 179)18111831) do not encode the RPI relations between coUinear and ultrasoft fields which 
can be determined using Eq. (fT6|) . 

The above results can be used to expand the RPI basis of operators in terms of operators 
in the homogeneous basis as in Eq. (H5l) . 



IV. APPLICATIONS 
A. Scalar Current 

As a first example to show how the expansion of a RPI current works, we expand the 
scalar chiral-even bilinear currents (LL+RR), for processes with a hard external vector q'^ 
up to order A^. In the basis built from superfields there is only one current that satisfies 
these conditions 

'^n,^J • (84) 
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All the other possible currents (for example ^ n,ui1fi(luGn''^^'^ n,d}2) have expansions that start 
at 0{X^) or beyond. To recover the basis with a homogeneous power counting, all we have 
to do is to expand (IMl) using Eq. (ISTj) . 

1 _ ^ 

'^n,d;J'^n,6j2 = ^^^Xn,LOi2Xn,uJ2 (85) 
O 7 ^Xn,OJa'^ Y-^t^i-^aT n-i-Xn,LJ2 7{ 7 r^Xn,u)iU-L'T''^T-i-'^a-i^2Xn,uja 

2uji[n-q)'^ 2uj2[n-qY 

1 _ 1 _ 

+ 7 77 — ^Xn,u^i^{q±-gB±) 

ljJa—L02Xn,LJa 

{uJi-uja){n-qY {uj2-uJa)[n-qY 

Xn.uj-i' ± ' (l±V^Xn, u)2 

1 — — Xn,w^1l>'yL-qi.Xn 



[n-qY uijj\ \n-qY OijJ2 

Thus all the (9(A^) terms (the twist-3 terms on the last three lines) are connected. Eq. 
agrees with the original derivation of these constraints given in Eqs. (122-126) of Ref. [16| 
The ease at which Eq. ( !85|) was derived demonstrates the power of the invariant operator 
formalism. In this example there is only one supercurrent to 0[}?\ so all Wilson coefficients 
are connected to the coefficient of the leading operator Xn,wi^Xn,<^2- Note that here all of 
the connected operators involve a g_L, which we have counted as (9(A°). We will see below 
that for situations with two coUinear directions, where in the end its natural to specialize to 
a frame where g± = 0, the connections tend to appear at higher twist. For situations with 
three or more collinear directions RPI will provide useful constraints on the basis already at 
lowest order. 



B. General Quark and Gluon Operators 

In this section we enumerate an operator basis for the general set of collinear quark and 
gluon operators up to 0{\^\ This basis is useful for many applications, and we keep our 
notation as general as possible. In particular we consider up to 4 distinct collinear directions 
(which for example could be used for e^e~ 4jets, or gg, qg, qq 2 jets). We also discuss 
a basis both for the homogeneous operators with a definite power counting, and for the RPI 
operators. 

For processes with a hard g^, the most general basis of homogeneous quark operators in 
SCET up to 0{X^) is 

Q{Oa) ^ - p Q{la) - p -pfa /n^N 

^ Xn-i_,iAj\ ^ a'Pn2± Xn2,uJ2 i Tq, (i^'^risx) Xn2,aJ2 ' 

Q{2a) ^ - p nV^",V^ , Y O^^'') = y F Y 

Ani,u)i^ afi ' ni±' n2± A.n2,uJ2 i ^ Ani,<^i-'- ap ' ni±' ni± An2,'^2 ' 

_ ^ r 10" IC^ -V n(2Q!) _ ^ -P -pfa ( „t2P \ 

^ Ani,!^!-'- a/3 ' n2± ' n2±'''-"2,<^2 ? ^ Ani,(^i-'- a/3 ' „-^_L ^^y»-'„3_Ly 



X 



n2,uJ2 } 



~ Xni,ajiF/3a (5'i3^3_|_)^^'P„2 Xra2,<^2 ? ~ Xni Fq,/? [7^„g {d^n-^^ ^,^1x1^2,(^2 5 

— Xni,uji'^al3 (5'^713±)lj3 (fi'^n4±)^^ Xn2,t<J2 ' ^ ~ {Xni,uji^lXn2,ui2^ {Xns, 013^ 2X714,1^4} ■ 



If we need to specify the subscripts we write for example 0^'^^\uji,uj3,uj4,uj2), with the uji 
listed from left to right. Due to the equations of motion in Eqs. (1471481) we did not need 



23 



to consider in ■ dnXn or gn ■ Bn- For each operator there may be a set of different Dirac, 
flavor, and color structures T^^...^,^ which depend on the particular phenomena being studied 
(including also two choices for color for the Fj in the four-quark operators O*^^*^). In general 
for each independent F^^...^,^ structure the operator has a Wilson coefficient that must be 
determined order by order in perturbation theory. We included in Eq. (1861) the mixed quark 
and gluon operators. For pure gluon operators up (9(A^) we have the homogeneous basis 

(87) 



(9 (Ob) 


= B^f" 




(9(l'i) 




"2 ,'^2 ' 




= B-^f" 


' n2± 112, Ld2 ' 


0(1/) 


= B^^ 


712 ,1^2 "3, 1^3 ' 


(9(2i) 


= B^^ 


-pfa pt/3 r^±u 
ni± ni± n2,LJ2 ' 


0(2.-) 


= 


ni± n2-L n2,uJ2 ' 


Q{2k) 


= B^^ 


' n2±' n2-L n2,oJ2 ' 


(9(20 


L' "i-L 


1 ^±u g±r 
'^?ii,a;ij '^ra2,'*^2 713,1^3 


Q{2m) 


= B^^ 


L' n2±'~^n2,u)2] n3,LU3 ) 




= B^'' 


g±u rpa 1 
"2,t^2 L n3±'^n3,aJ3j 




= B^^ 


f2-^i^ f2-^^ f^-^^ 









Here we do not need to consider operators with gn-Bn and gn-dnB'^j_ because using the 
equations of motion in Eq. fHHl) they can be written in terms of the operators in Eq. fISTI) . 
and are hence redundant. 

To setup the computation of constraints on Wilson coefficients we also need to build an 
RPI basis of operators using the objects in Eq. (129|) and id!^. Because each operator will 
be RPI, its Wilson coefficient is truly independent of those for other operators in the basis. 
The RPI operators can then be expanded in terms of homogeneous operators made out of 
of gauge invariant objects, and doing so we obtain operators in the homogeneous basis with 
all the constraints coming from reparametrization invariance. The number of constraints on 
Wilson coefficients is equal to the number of homogeneous operators minus the number of 



RPI operators, once we have accounted for linear dependencies |35!, [36 . 

Let's construct the RPI basis of operators which is the analog of those in Eqs. (!86l) and 
(!87|) . The operators with no id^ derivatives are 

where the basis of Dirac structures F, and contraction of indices fivar in Q(o») depends on 
the kind of current we are studying. For cases without a the subscripts uji are erased 
and RPI operators are multiplied by the Aij factors shown in Eq. (!33|) . Recall that we 
do not have a good power counting in the RPI basis, this basis makes the RPI properties 
transparent but the power counting more tricky. When Q^"-'^^ and Q^'^^) are expanded in 
terms of operators that are homogeneous in the power counting, they contain a leading 
order term, so they are relevant operators to consider at LO. Of the RPI objects only id^ 
starts at leading order, so theoretically we can construct an infinite set of LO operators using 
{idn)^ for any k. However, the structure of this operator provides additional constraints. 
In particular the 0{\^) term is id!^ = {n^/2)Vn + and the collinear momentum P„ 
acting on a n-coUinear field such as just gives a number, lji, which can be absorbed 
into the Wilson coefficient C{ui,U2)- For cases with a this implies that adding id^'s 
in a scalar operator (where all vector indices are contracted) most often gives an operator 
that differs from one we already have only at 0{X). For these scalar operators we can count 
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id^ ~ 0{X) when determining which RPI operators are required, and for simphcity we follow 
this counting in the remainder of this section. If we have an operator with a free vector 
index fj,, then this index can be carried by id!^ = {n^/2)Vn + ■ ■ ■ , and the partial derivative 
does count as 0{X^). 

The expansion of the RPI operators in Eq. (I88l) in terms of homogeneous operators up 
to C(A^) is 

Q(Og) ^ ^(1) ,^{2) I r^(2)^ (gQ\ 

^ ni.oji n2,ci;2 "i,<^i "2i'^2 "i,<^i "■2i'^2 V / 

' nijUJi n2,uj2 "li'^l 1^2,1^2 n2,uj2 V /' 

(JT 



(Og) ^ QW^^v q(1)ctt q(2)iiv q{1)ut ^(I)mi' q(2) 

"ni,a;i "n2,a)2 J^liiil "n2,(i2 "n2,il'2 



where the and are given in Eqs. fIST]) and fl5^ . To look for RPI relations the results 
of this expansion must be compared to power suppressed operators which also can generate 
0{\^) and (9(A^) terms. Up to this order the power suppressed operators involving two or 
more quark fields are 

Q{ia) ^ xfy - r 7(9" v|> . O^^''^ = vi/ - r id"" v|> . fqn'i 

q(1c) ^ ^ . p (-/3/3' ^ _ Q(2a) ^ - F / ^ - 

Q(2fc) ^ ^ - r / z^" \l/ - O^^''^ = ^ - F 'id" id"' ^ - 

q{2/) ^ ^ . p C'^^'- ^ - O^^^^ = ^ - F aC"^- C'^'^ - ^ - 

^ ^"iii^i-^ app ='n3,ii3 '^'^n2 ^"2,1^2 5 ^ ^ 0/370 ^n3 ,1^3 ^714,1:^4 ^ 't-2,'^2 5 

q(2/i) = - Fi\l/ - 1 - F<,\1/ - 1 

Again a minimal basis for Dirac structures F will depend on the process being studied 
and may differ between the various Q^*^) operators. Such a basis will also in general differ 
from the one for the homogeneous operators in Eq. (!86|) . We will adopt notation such as 
Q^^^''(a)i, cjs, UJ4, UJ2) when we wish to specify these subscripts. For a field basis for the higher 
order operators with gluon fields (whose expansion starts at 0{X^) or 0{\^)) we have 

Q{ld) Qfiu ■r.a Qar Q{le) QfJ.u ^a q<tt 

^ ^ni,Cji''^n2^n2,uJ2 ■> ^ ^ni,(ii '^ni yn2,a;2 ' \^^) 



q(1/) ^ nt^'^ G'^r Qaf3 Q(2i) Q^lu -ga • o/3 nc 

^ ^ni,(ii^n2,a;2^n3,(i3 5 ^ ^ni,(ii ''^n2 ''^n2 ^r!_, _ 

^ 5='ni,wi'''^n/<^ni 5='n2,il'2 ' ^ 5='ni,iii '''^ni '^^n2 5='n2,tJ2 

Q(20 = ^M^', l^<xr^ ^a/3_ Q(2m) ^ ^^^^ _ ^-gy gar ^ 



,1^3 ' 

Q{2n) ^ ^/ii^ ^o-T_ r-o7 1 QC^"^ = G^" - G"^ G'^^ 

We will include a basis of Dirac structures and expand the RPI operators in Eqs. (MJ|) 
and flOTl) in terms of the homogeneous ones in several of the examples below, and consider 
whether there are non-trivial RPI relations on a case-by-case basis. 
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C. Deep Inelastic Scattering for Quarks at Twist-4 



In this section we consider spin-averaged DIS at twist-4. This provides a test of our 
technique of constructing a minimal basis, for an example where the basis is already well 
known 20|, l2l|, |23|. We will see that RPI constrains the Wilson coefficients of the homo- 



geneous coUinear operators. Our analysis is really of scalar operators with one collinear 
direction, g_L = 0, with overall derivatives set to zero. DIS is the most popular application 
for these operators, so we frame our discussion in that language. For simplicity we consider 
the QCD electromagnetic current = q'j'^q for one-flavor of quark. (We briefly discuss 
the generalization to non-singlet operators in a footnote.) The study of higher twist in DIS 



and related processes is an active field of research, for example 37|, |38|, |39|, |40|, |4l|, |42] . In 
the language of SCET, DIS was first studied in [3], whose notation we follow. The virtual 
photon has momentum transfer = —Q^, and x = Q"^/ {2p ■ q) is the Bjorken variable. 

In the Breit frame the momentum of the virtual photon is q'^ = Q{n^ — n^)/2, and 
the incoming proton momentum is = n^n ■ p/2 + n^^m'^/{2n ■ p) where irtp is the mass 
of the proton. Expanding in ntp/Q we have n ■ p = Q/x — xrn^/Q + . . .. The energetic 
proton has a small invariant mass p^ = nip ~ ^qcd' ^^"^ Breit frame it is described by 

collinear fields in the effective theory with a power counting in A = Aqcd/Q- It is convenient 
to pick this frame in order to be able to assign definite power counting to momentum 
components. What reparametrization invariance enforces is that all results are invariant to 
small perturbations about this frame, encoded by changes to the collinear reference vector 
n^. Since these changes are small we are free to use the same power counting when studying 



the RPI relations. There is a larger class of frame independence, which says for example that 
the same results would be found if we compare an analysis in the Breit-frame with an analysis 
made about the initial proton rest frame, but this set of "big" frame transformations does 
not encode non-trivial dynamic information that relates coefficients of operators at higher 
twist. All final results are of course entirely frame independent. 
For spin-averaged DIS the hadronic tensor has the structure 



T, 



+ '-f) Ux,Q^) + {p, + I) {p. + I) T,(x,Q^), (92) 



where 



T,Ap,q) = lj2^P\T,.{q)\p), f^u{q)=i J rfSe^^-^T[J,(^), J.(0)] . 

spin 



(93) 



The scalar structure functions Tj can be projected out of T^^, using 



T,{Q\x) = --[g 



4x^ 

Q^ + Amlx^ ^''P'j^M-' 



2x2 



12x2 



+ ArripX^ 



(94) 



The expansion of Ti and T2 has been carried out up to twist-4 with the Wilson coefficients 
determined at tree level in Refs. 20|, l2l|, |23|]. To simplify our calculations we will make use 
of the fact that the projections in Eq. (jUj) commute with taking the proton matrix element. 



26 




Here [duJk\ = duoi ■ ■ ■ duo^ is the integration measure over the independent parton momenta 
ujk carried by the Wilson coefficients C]*^ and the operators Oj. The superscript [i\ indicates 
that the Wilson coefficients for the two tensor structures will in general differ. We also 
consider a basis of RPI operators Q^- by writing 



Unlike the Oj's the Q,'s do not contain contributions of a definite order in the power 



coefficients of the Oj's. A connection would mean, for example, that the one- loop 
coefficient for a twist-4 operator is determined by a coefficient at twist-2 at all orders in a^. 

We first write down a gauge invariant basis of chiral-even quark operators that are ho- 
mogeneous in the power counting. This can be done using the general basis in Eq. (|86|) 
with all directions = n. Furthermore, since the DIS matrix element is forward, we have 
{p\\P^O]\p) = for any operator O. Thus we are free to integrate ±-label momentum 
operators by parts, and hence can ignore all terms with V\_^s in Eq. ( 1861) . (If we consider 
our analysis to be of the general scalar operators with one coUinear direction, then this is 
the only simplification that we make which relies on the form of the final matrix element.) 
For simplicity we also drop the square-brackets from inside O*^^-^) in Eq. (15^ . A minimal 
basis of chiral-even parity-even Dirac structures between the ?7,-collinear quark fields is easily 
constructed using the properties of the SCET Xn fields. We have i) just {^} when there are 
no vector indices on fields, ii) no elements at all when there is one vector index, and iii) just 
{'^g'^ , or {i^g^ , "^7x71} foi' two vector indices on fields. Here ii) is the standard 

fact that the spin-averaged case does not have twist-3 terms. (For polarized DIS it does 
not suffice to only consider the scalar operators.) For the four-quark operators we can have 




3 



(96) 
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Fi ® r2 = ® '^75 ® "^75} and color structures 1 (g) 1 or T^. Thus the basis is 



^1 Xn,(^i 2 Xn,uj2 } 

n - t 

^3a Xn,uji 2 



Xn,iL}\ 2 Xn,<^2 



(97) 



04a = Xn,Loi - {g^n±)^aT±Xn,uJ2 ; 



09 = 



Xn,uj\ 2 Xn,aJ2 
Xn,uii 



^7.T^Y 

'1 2 'S-' An,a;2 



Oib = Xn,u;i-T±i9^n±)^3Xn,uJ2 ; 

= Xn,c.i 2 Tr [(^^„^)<^3(^^„^)a;4]Xn,a;2 > 
— Xn,uji Tr [(5'^„_|_)aj3 (fl'^/x )t<J4] Xra,a;2 ' 

r 1 r_ ^ 1 

Xn,aJ3 ^75Xn,aJ4 



^8 /\n,uji 2 

Xn,uii ■^75Xn,a;2 
A.n,uJi 2 An,ci;2 



Recall that in an operator like O2 the position space analog of is to translate all gluon 
and quark fields in Xn,a;2 in differentiate twice with respect to x^, and then set x_l = 0. 
The basis shown in Eq. fl97|) can be used to describe twist-4 effects in DIS at any order in a^. 
Note that we have already discussed and taken into account the quark and gluon equations 
of motion in the general result in Eq. fISBl) and hence already in Eq. For O5J there 

are two color structures associated with the product of i3„'s, but these are picked out by 
consider Wilson coefficients Csj that are odd or even in the exchange cjs ^ 074. The forward 
proton matrix element of these operators will be proportional to an overall (5-function, which 
is 5{uji - LJ2) for Oi,2, 5{u}i + c^s - ^2) for 03a,3fe,4a,46, + c^s + u;4 - U2) for Os-s, and 
5{ui + Lo-s - - ^i) for 09-12- 

Next we derive the analogous results for the RPI basis of chiral-even operators. From 
Eq. (jnSD the hadronic tensor operator T^^ depends on g'^ which we use as our reference 



vector. To construct this basis we cannot use or n^'. Comparing Eqs. fl^?5]) and Eq. flM|) 
we see that it suffices to construct a basis of scalar operators for the expansion of g^^T^y 
and p^p^T^y. The forward proton matrix element of the expansion of these operators then 
yields an expansion for the observables Ti and T2. Thus, for the scalar basis we allow any 
number of g's to appear, but only zero or two p's. This implies that at twist-2 there is only 
one RPI bilinear quark operator 



Ql = ^n,<ij^n,c: 



(98) 



At twist-3 there are no scalar chiral-even RPI operators. The candidate operators \l/n^^n^n 
idn)'^n„u2 ruled out by the equations of motion in Eqs. fl^ and (15^ . 



and ^n„c^i„_ _ „ . . . , ^ 

Another possible operator is '^n'pip ■ dn)'^n but it starts at twist-6, since (p ■ 9„) ~ O(A^), 
being suppressed either by an n-p or n-(9„, and |^ adds another factor 2 to the power counting 
when it is squeezed between the n-collinear fermion fields Xn- All the operators with Q!^^ ., like 
for example ^nlfiQuGn'^'^n, have expansions whose lowest term is twist-4 because the Dirac 
structure of the twist-3 component of this operator vanishes between the n-collinear fermion 
fields, since Xn'4Xn = 0. Thus the power suppressed terms start at twist-4 in agreement with 
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the homogeneous basis in Eq. (197|) . Writing out the RPI operators different from zero at 
twist-4 and not connected by operator relations we have 



Q4 = -9^^n,^Jl^,quG'^%^laqpGn^4^n,c;J2 , Qs = -9^^n,^Jlf.qu1aqpTT[g!;^%^g^'^^J'^n,^2 , 

Qe = -g^^n,^JqMGn,^3TaGZ,^n,^2 , Qt = -^?'^n,ii^ ^M^-Tr [(6;„,^3 )'^„e?°'^J *„,^, , 
Qs = [^n,^r0n,^2] [^n,u:i0n,^,] , Qg = nMH5'^ n,^2][^ n,^Jl5'^ n^] ' 

One can think of other possible operators at twist-4, but all of them are either ruled out 
by the equations of motion and operator relations, or start at higher twist. For example, 
there are not operators with both and Qj^'^, like ^'n,ii7^jPi/^^'^3^n,cJ2 5 because they all 
start at higher twist. We have integrated by parts making all derivatives act to the right, 
since here our interest is in forward matrix elements, and we removed id^Q^cj with the gluon 
equation of motion in Eq. ([S3D- The operator ^n,wAi.^<^n - idn)'^n,(b2 = ^n,ij^^n#n^n,i2> 
and is removed by the quark equation of motion in Eq. (ISTl) . For the operators with two Q's 
only the structures in Q4_7 have expansions that start at twist-4. For example, Q'^%^Gn,c:j4,au 
at LO is proportional to {gB^)^^{gBj^)^^n^^vy so closing the indexes with 7^ or 7^^ generates 
a ^ that next to Xn gives zero. 

It is less obvious that operators with one Qn and one idn are redundant and can be 
eliminated from the RPI basis. Consider the operator 

Q* = ^n,C.Aq^^^9Qn,C^,^d:-^n,C.2 ■ (100) 



To remove it we use a manipulation discussed by Jaffe and Soldate in Ref. [2l|]. First we 
write 

g^'^td: = GH'^tv: - ^r[i/(? ■ ^dn)qat9g:.] , (loi) 

and note that the term with two ^„'s can be ignored since it is already in our basis. Next 
using the definition (!36l) we can write 

q^igQil'tV: = q^[zV;,zV:]tV: = l{-q^[iV:, liV;,zV:]] - {zb^)Hq.dn + iq-dn{ib^f] . 

(102) 

The double commutator term is turned into a four-quark operator by the gluon equations 
of motion in Eq. (l53l) . For the remaining terms we write {iDnY = iTpniTpn + f 0"/^^^^'', where 
the a^y term gives Qg, and terms involving {i'pnY are turned into the operators Q2, Q3, 
Q4, and Qg by the quark equation of motion in Eq. (1^ . (They are not simply set to zero, 
since [iTprl\ does not commute with 5{uj — 2q ■ idn).) Finally, we can also rule out the only 
other non-trivial operator ^ n,LsAipn1 iJ.qyQ^L33^ n,L32- Using the gluon equation of motion we 
write 

^n,u^Ml,quGZ~,'^n,u^2 + ^nMl.^^GZ^.qu'^n,^^ = -2Q* + • • • > (103) 
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where the elhpsis denotes operators with two QnS or four-quark fields that are part of the 
basis. The Bianchi identity in Eq. (|54l) gives another relation for the two operators on the 
LHS of Eq. (11031) and implies that they can be written in terms of Q3, Q^, Q4, and Qg. 
Thus both the operators ^n,,:uJiq^nr,iqiyQn%g'^n,^2 and ^n,&JltiiPnQn%^qu'^n,^2 are redundant. 
Finally we note that the order of the QnS in an operator like Qg is not important, since 
we can always symmetrize or antisymmetrize its Wilson coefficient in 1^3 and 1^4. Note that 
when considering the transformation of the operators under charge conjugation one must 
consider both the operator and its Wilson coefficient. We discuss an example below in 
Eq. ([no]). 

The number of independent RPI operators in Eq. (|971) is smaller than in the basis of 
homogeneous operators in Eq. ( |99l) . implying that there exist further constraints on the 
Wilson coefficients of the homogeneous basis at twist-4. To find the constraints we must 
expand the operators in Eq. fl99|) in terms of those in Eq. fl971) . We start with Q4 through 
which are in one-to-one correspondence with operators in the homogeneous basis, 



Q4 
Q7 

Qio 



a;3a;4 
4{n-q) 

-UJ3UJ4. 
A{n-q) 
1 



O5, 



Q5 

Qs 
Qii 



A{n-q) 
1 

[n-qY 
1 

{n-qY 



On. 



Q9 



ci;3u;4 



Or- 



A{n-q) 



(n-g) 



(104) 



Here the order of the Ui subscripts in operators on the left exactly matches up with the Ui 
subscripts on the right. For the remaining operators whose expansions start at twist-4 and 
for Qi that starts at twist-2, we have 



Qi(a;i,c2;2) 



{n-qy 



{n-qy 



^Oiai^^l, ^3, ^2) + ^04b{uJl,UJ3, UJ2) 
- OsaiuJl, UJ3, UJ2) + OsbiuJl, UJ3, UJ2) 

— Oiai^^l, ^Z, ^2) + O^bij^i, t^3, UJ2) 
Ll)2 



(105) 



n-q 



+ 
+ 
+ 
+ 



Oi{uJi,uj2) + ^ <! ■ 

n-q [n-qY [ L UJ1UJ2 

2 d 2 



r -1 



+ 



d 1 



+ 



+ ... , 

d 1 



+ 



02{UJI,UJ2) 



{uJa-UJ2Y duj2^^a-^2 

-2 d 2 



duOi UOi dU2 002 ■ 



(u;i-U;a)2 duJxOJx-OJa 

Oia\yJa,0Jx-UJa,0J2) + 



ijj\ijj2 dijj\ 
-1 d \ 

U)\U)2 duJ2 OJ2 



03a(cUa,t^l-t^a,t^2) " O^bioJ a,OJx-UJ a, OJ2) 

d 1 



duJ2 iOa 

O4b{0Jl,UJa-UJ2,UJa) + 3 O^bi^Ja, OJl-UJa, UJ2) 

dUi LOa 



+ 



Here the ellipses indicate terms involving operators 0^-i2 that have already occurred in 
Q4_ii and hence they are no longer important for determining the linear independent combi- 
nations. It is interesting to note that expanding the operator gives the same combination 
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of Osa and Oaf, that appears in Q2 — 4c<Ji/co'3Q3, so even if we had not ehminated from 
the RPI basis, the imphcations for the homogeneous basis would be the same. 

The three RPI operators in Eq. (11051) have expansions in terms of six homogeneous 
operators Oi, O2, O^a, O^b, O^a, and 046, so there are three RPI relations. The Wilson 
coefficients of these six homogeneous operators are determined by three coefficients, Oi,2,3 
in the RPI basis. It is convenient to trade C*i,2,3 for the three coefficients Ci, C^a, and C^i,. 
The remaining coefficients O2, C^a, and C^b are then determined by RPI. We find 

n-q r 1 Id Id 



02(^^1, = < \ 3 \ 3— \Ci{uJi,UJ2) , 

n-q IUJ1UJ2 uji duji uj2duj2J 

Ciai!^l,^Z,^2) = -7;C3a{i^l,i^3,i^2) " 7^ Csb{uJl, UJ3, ^2) H Ci{uJi, UJ2-UJ3) 

2 2uj2 n-quj2UJ3 

n-q{uJ2+uJs) 

Ci{uJi+UJ3,UJ2> , 



n-quj3{uj2Y 

n ( \ ^'^ n ( \ "^n ( ^ i ^"^l^i-^s)^. x 

L^bK^l,U3,U2) = -7; C3a{i^l,UJ3,i^2) " 7;'-'3b[uJl, UJ3, ^^2) H 7 -^p Ci[uJi, UJ2-UJ3) 

I n-quj3[ujiY 
^ 01(^1+^3,^2). (106) 



n-quJiUJ3 



We have cross-checked the relation for O2 with a tree-level matching computation. Note 
that 02(0^1,072) multiplies a matrix element that gives 5{uji — 0J2), while O4a,4b(oji, c^s, o;2) 
multiplies a 5(ti;i + 0^3 —0^2), and that we have used these 5-functions at various intermediate 
steps. That is, the result in Eq. fll06p applies for a basis of operators, whose matrix elements 
have vanishing total derivatives. 

Our operator bases can be compared to the flavor singlet and parity even basis of Jaffe 
and Soldate in Ref. 2l|] which has one operator at twist-2, and 12 operators at twist-4.^ 



There is a simple correspondence between the 11 operators in our RPI basis in Eq. fl98|99p 
and the QCD operators in their basis. The correspondence is one-to-one for Q^, the four- 
quark operators Qg-n, and the operators Q2 3 that have one Qf;^^ . For the operators with 
two Ql^'^^s we have four operators compared to their six, but the difference is accounted for 
by the way in which the twist towers are enumerated. We used continuous a>j's where even 
and odd symmetry under the interchange UJ3 ^ 004^ encodes two possible color structures 



with f^^'^ and d^^^ , while Ref. 2l| uses a discrete basis with integer powers of {in ■ Dn 



where the choice of which operators to eliminate by integration by parts implies that the 



^ The notation in Eq. ([97]) suggests that all quark bilinears are flavor singlet contractions if Xn has multiple 
flavor components. To incorporate other possibilities for the flavor indices is straightforward (2l|. We 
consider Xn as a doublet of SU(2) flavor, or a triplet of SU(3) flavor, with elements Xn- For photon 
currents one has a charge matrix in flavor space in each QCD current, which is Q = diag(2/3, —1/3, —1/3) 
for SU(3). Thus, at leading order in the electromagnetic interactions one must simply introduce a in 
all bilinear-quark operators, through in Eq. ([57)1 . When counting the four-quark operators to 
O^^ induced by photons we double the number of operators because there are two possibilities, ^ 1 
and Q (S) Q- In this notation the flavor singlet contraction for the four-quark operators is 1 (g) 1. For the 
RPI basis of operators the analysis of flavor structures is identical, and hence flavor does not modify the 
constraints in Eq. p06p . 
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two color structures yield different operators. Our homogeneous basis has 14 operators up 
to twist-4, and most closely corresponds to an enumeration of an operator basis in terms 
of the so-called "good" quark and gluon fields. The good quark and gluon fields have been 



discussed in Refs. 2J, |43|, |4J] . In this basis the power counting is manifest. From the three 



RPI relations in Eq. (I106P the number of independent short distance Wilson coefficients 



is 11, and so encodes the same amount of information as the OPE basis from Ref. [21 
Note that there is no room in the traditional OPE in DIS for a correspondence with higher 
order operators with ultrasoft fields. In our language, the validity of the OPE for DIS with 
generic x implies that ultrasoft degrees of freedom are not needed, and one can consider 
that fluctuations from that region are reabsorbed into the coUinear fields. 

When the basis of bilinear quark operators is considered in the forward proton matrix 



element it can be reduced even further as discussed in detail in Ref. [23|. In this process 
it is found that the matrix elements of operators like O2, O^a, and O^b do not provide 
independent information. Hence at this level the RPI relations in Eq. fll06p do not appear 
to have practical implications. 



D. Deep Inelastic Scattering for Gluons at Twist-4 



Next let us consider the minimal basis for pure gluon DIS operators up to twist-4. We 
proceed in a similar manner to our construction for quarks, first writing the homogeneous 
basis and then the RPI basis to check if reparametrization invariance provides constraints 
on the homogeneous operators. The homogeneous basis is 

(107) 





= Tr 
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U{gBl^UVl{gBi^ 




05,6 


= Tr 


\9B'n± 
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07,8 


= Tr 




U9B:^U]TT[{gB:^ 
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= Tr 









where F 



1,2 

fiua[3 



{9fiu9a/3, gfiaguis} and the traces are over color. Recall that the equations of 
motion (l48ll were used to eliminate the operators gn ■ Bn and in-dn{gB'^^). Again since the 
basis is designed for taking forward matrix elements we are free to integrate by parts and 
hence we do not consider . There is a third tensor structure, F^^^^ = g^pgau, that can 
also be considered for 0^-%, but which can always be eliminated. For 03^4 this is done using 
integration by parts and the cyclic trace, giving 



'^[{9B:^)^A9BW).,Vl{gBij^,]V^ 



-Oi{u2, ^3, ^^l) - 03(u;3, tUi, U}2) 



(108) 



For 05_8 the cyclic property of the trace suffices to eliminate Fj^^^^ in an analogous manner. 
The operator Ti[{gB':^^)^-^V1^[^gB^_^)^^{gB^j^^^ is also not needed in the basis because it 
can be put into the form of the operators O3 and O4. This is done by acting with the 
on the two B±_'s to the right, using the cyclic property of the trace, and again noting that 
03^4 encode all orderings for the Ui subscripts. 



32 



For forward spin averaged matrix elements the RPI basis of gluon operators up to twist-4 

is 



Qi = gi.g/3Tr [igGil^^igGn'^^^] g^,a , (109) 
Q2 = Tr [igQn^^^igQn'!u^ g^^agu/3 , 

Q4 = QuqpqaTT[igg^'^^Jgg"^^Jd^igg^'^^^]g^p , 
Q5,6 = ?<^?/39agATr[i^^^':^i^^°J^ic/^^"^^i^^;^Jr^'2^^ , 
Q7,8 = ?<^?/3g-?ATr[i^^r^i^^^JjTr [i^^^"^3i^^;^J^Jr^'2^^ . 

Here we remove a possible operator gi^g^gTr [(i(7^„^ij)'^(i9)^i(y'^"^^] by writing {id^g'^'lj^ = 
iidfj_)id^g"^^^, and then using the Bianchi identity in Eq. ( !54l) to rewrite this operator in 
terms of operator with two ^n's, plus (^<9^)*c^"^f^2 and {'id^)id^g'^"^^. The last two terms 
are removed by the gluon equation of motion. There is no need to include the analog of 
Q4 with the idf^ acting on iggnlj^^ because it is related to Q4 by integration by parts up 
to a term, id^j.igg'^^^ that reduces to other operators through the gluon equation of motion. 
Again the cyclic nature of the trace allows one to remove T^^apr Qs-s- 

In order to consider the effect of charge conjugation on these basis one must consider the 
transformation of 

'|dci,]C.(^,)Q.(i,). or /|d.ja(.,)0.(.,,). (110) 

where Ci is the Wilson coefficient associated with Qj, and Cj the Wilson coeffi- 
cient associated with Oj. We can impose constraints on Ci{ujj) and Ci{ujj) such that 
(11 101) is C-invariant. For example, note that under charge conjugation Q3 transforms 
into -gi/gaTrfisf^^lllgiisi^^^^^iiSi^^^^Jsi^aSi^p, so to make it C-invariant we impose that 

(73(0)1,0)2,^1)3) = — C*3(o)3, 0)2, 'i'l)- Similar considerations apply to the homogeneous ba- 
sis. For example, the combinations 03(o;i, 0^2, cf3) — 03(o;2, oJi, 0^3) and 04(o;i, 0^2, ti;3) + 
Oi{ui, 003,002) + 03(0^3, 0^2, cui) are even under charge conjugation. 

Next we must expand the RPI basis in Eq. (11091) in terms of the homogeneous basis in 
Eq. (11071) to find possible constraints. We first expand Qs.g, they have only operators with 
four gB'^'s, that is Os.g, 

L01LO2LO3UJ4: LO1UJ2UJ3UJ4 
^5,6- ^ ^5,6, ^7,8- ^ <^7,8- iJ-iij 



Next we expand Q3 4 to find 

Qs = A^^^'^ \ \ ~ Oa{uji, t02, UJ3) - 04,{lo3, LOi, LO2) - 03(o;2, o;3, o;i)l H , 

4(n ■ q) ^ 

Q4 = ^^^^^^ 04{uJi,uJ2,uJ3) - — 03(o;2,o;3,o;i) +..., (112) 

o L UJi J 

where we integrate over the repeated Ua variable. The ellipses in Eq. (I112p indicate terms 
involving operators Os-g that have already occurred in Q,^_g and hence are no longer im- 
portant for determining the linear independent combinations. Eq. (11121) implies that O3 and 
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O4 have Wilson coefficients that are independent of other operators in the basis. When we 
expand the remaining RPI operators 2? we may also have terms with Oi,2,9 which have 
two gB'^^'s. We find 

Q,(-„ = {u^^u,) + ^ (2 - — - — .2 j O (u;,, .2) + . . . , 

Q2(t2'l,CJ2) = . ■ . , (113) 

where the ellipsis indicates terms involving operators O^s that have already occurred in 
Q3_g. The fact that does not occur in the expansion of any of the RPI operators indicates 
that it is ruled out by RPI (explaining why we listed it last in the basis). Furthermore, the 
operators and only enter in the combination obtained from expanding Q^, and so 
their Wilson coefficients are related by 

C2{UJI,UJ2) = [2 + UJi- \-U2- 



n-q\ duoi du2^ UJ1LJ2 

n-q f d d 



(^^A7 ^^^^)'^^*^^^'^^'* ■ ^^^^^ 



uJiUJ2n-q\ duJi 

For the gluon DIS operators the RPI relations are similar to that for the quark basis, 
namely it is the collinear operators with P^'s that are constrained. This was also observed 



in Ref. for the heavy-to-light currents at second order in the power counting. Overall 
there are eight homogeneous operators for spin-averaged gluon DIS up to twist-4, and seven 
independent Wilson coefficients. 



An analysis of twist-4 gluon matrix elements was done in Ref. 451] using leading-order 



Feynman diagram, based on the methods of Ref. [23|. To the best of our knowledge, the 
complete linear independent bases of twist-4 pure glue operators given in Eq. fll07p and 
(11091) have not been given earlier in the literature. 



E. Two Jet production: n-n' operators 



An important application for operators with two-coUinear directions, n-n', is the study 
of two jet phenomena and event shapes. The effective theory SCET has been used to 
study jets at leading order in the power expansion and various orders in the expansion 



in Refs. 27, 46, 47, 48, 49, 50, 51, 52 
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54J . |55| . Another interesting application is to 



describing parton showers with SCET 56|, |57|, where both leading and subleading operators 



with two-coUinear directions play some role. In this section we study the leading and first 
power suppressed quark operators with two-coUinear directions. For two jet processes it is 
convenient to use the center-of-momentum (CM) frame where the two jets are back to back. 
In this frame we can take n' = n so that n'-n = 2. Our main interest will be in the operators 
that do not vanish in this frame, however part of our discussion touches on the additional 
operators that do. 

To be concrete we consider operators that appear in two jet production from a virtual 
photon of momentum q^ in e"^e^ JnJn'- In QCD the fundamental hadronic operator is 
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the current = ip^f^ij;, which is conserved d^J'^ = or q^J'^ — 0, is odd under charge- 
conjugation, and transforms as a vector under parity and time-reversal. To describe high- 



energy jet production this current is matched onto a series of SCET currents J] 
with Wilson coefficients Ce{uJi), 



oo „ 



2-jet 



(115) 



Here k denotes the power in A, the subscript £ denotes members of the basis at a given order, 
and the uji are the set of gauge invariant momentum fractions upon which the operator 
depends. We also sum over all coUinear directions n and n', and the appropriate ones for a 
given computation are picked out by the jet-momenta in the states. Because of this sum we 
are free to swap n ^ n' when considering symmetry implications. The C, P, and T symmetry 
properties of Ce,{(jJi)jf'\uJi) are the same as J^, and they also satisfy current conservation, 
q^[jf^\uJi)Y = 0. Finally, since the matching takes place at a hard scale where perturbation 
theory is valid, the SCET operators should have the same LL + RR chirality as J^. 

We first construct a basis of SCET operators that is homogeneous in the power counting 
and with even chirality. For the construction of this basis it is convenient to define 



uu 

9t 



9 
n-q 



q^q" 



g2 
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■n!^ , 



— n'^ H n^" 



n -q 



where r1 is odd under n n' and r+ is even. We also define s'^ as with n ^ n and 
n' — > n' . Four of these objects are transverse to 5'', q,j,gj^ — 0, q^^^ — 0, and q-r- — q-S- — 0, 
which is helpful for satisfying current conservation. For constructing the homogeneous basis 
it suffices to consider the vectors {r_, g, s_, s+j in place of {n, n, n', n'}. When we specialize 
to the CM frame, g^_|_= ?^'j_ = 0, s± = Tr±, and the vector r1 = q^, and hence and s± do 
not need to be considered. 



In a general frame the LO operator is Xn,ujJ^^Xn,w2 with F'' = {7^, r^! ^, r^! /_ 



flV -I- 

I 9t 



r'tfjr,gi^r^f-,gi^r^f+} plus terms where r+ or r_ are replaced by s±. No terms with q^ 
are allowed by current conservation. Things become much simpler if we focus on operators 
that are non-zero in the CM frame. In the CM frame x 



0, 



x^ 



-Xn. 



■^2 



0, 



and the vectors r+ and s± become redundant, so there is only one operator at lowest order 



(117) 



Here cjj = {c<Ji,c<J2} and for brevity we suppress the index on the LHS. 

To construct a homogeneous basis at NLO we again consider only operators which are 
non- vanishing in the CM frame. In the CM frame we can take the total transverse momentum 
of the jet equal to zero, so we have the relations Xn\u}J^ tjiP lXn,uji = Xn' ,uij^ jxV^x^ Xn.uji = 0, 
with F^ any gamma structure, and hence do not need to consider operators with a single 
V^. Again all operators with a ^ or /_ vanish, as do those with q ■ {gBn±) and r_ ■ {gBn±), 



and the analogs with n 



n 



Operators with three 7's can all reduce to operators with a 
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single 7 plus terms that are zero in the CM frame. This implies that at NLO there are only 
two operators 



Linear combinations of these two SCET currents can both be made odd under charge con- 
jugation by imposing appropriate conditions on their coefficients under ui <-> —uj2- 

To see if there are constraints on the Wilson coefficients we write down a basis of RPI 
operators up to NLO. The objects 7^ and g!^'^ are invariant under RPI and can be used for 
this construction, but the object r'^ cannot. We find the basis 

= ^nVl7^^n,^2 , (119) 



n,uj2 ■ 



Here we do not write down RPI operators which vanish in the CM frame when expanded, 
such as ^ n' g^u'^d'^'^ n or operators with only the Dirac structure ^. This set also includes 
three 7 operators since in J^^"* replacing T^ap = gj^alis T^a/j = 7^7^7/3 gives an operator 
that vanishes in the CM frame, and any other order for the 7's is then redundant. The same 
is true for T^^fs = ijiilalp — 'f'gjialp- Analogous arguments rule out three 7 terms replacing 
the tensor in 3^2^ . There are no others operators with id^ or id^, besides Ja^g To see 

why, notice that for the operators with id^ and only the contraction with gj^y^ has the 
potential to give a LO term. Momentum conservation requires = id'^ + id^, and because 
q^gj^x = 0, we can exchange id^, and id^. The operators Jg^g correspond to keeping id^ 
when we have a ^"a, and id^, when we have a ^"f-,. 

The number of operators in Eq. flll9p is greater than that in the homogeneous basis, and 
when expanded Jo,i,2 Jo,i,2- Thus the operators in Eqs. (11171) and flllSp are not connected 
by RPI. For two jet production the constraints imposed by considering the CM frame are 
strong enough that RPI provides no further information. (RPI could still constrain the 
homogeneous basis of operators in a general frame, but does not have practical implications 
for determining the basis of operators for an analysis to be carried out with homogeneous 
operators in the CM frame.) 



F. Three Jet Production: ni-n2-n3 operators 

Here we analyze operators for three jet production. As in the two jet case, we consider 
production of jets from e~^e~ scattering through a virtual photon. To construct the minimal 
SCET basis needed for a matching we could proceed like in the previous cases, by writing 
down both the most general homogeneous basis and RPI basis consistent with the symmetry 
of the process, and expanding the RPI basis to find possible connections. In the two jet 
processes the interesting terms in the homogeneous basis are made of only two operators 
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up to NLO. However, for three jets the homogeneous basis has many operators at LO since 
we have three distinct directions rii, 77.2 and 713. With only two directions we could greatly 
reduced the number of operators by focusing on the ones that do not vanish in the center 
of momentum frame, meaning those that do not vanish when rii —>■ n, n2 —>■ n, and q± — 0. 
This choice rules out many operators because of the relations Xn^±i'Xn = Xn^±^Xn = 0, 
where r_|_ is a Dirac structure without t/l 01 iji factors. With three directions there is more 
freedom, for example the perpendicular direction of ni-fti is not the same as the one of 
n2-n2, or of n^-n^. Now to construct the homogeneous basis, we can use ni, ni, 77,2, n2, ^3, 
ns, 7^, so in the three jet case we have a bigger set of objects available. 

On the other hand the RPI basis still has a reasonable number of objects, namely \E'„-, 
(9(^., 7^^, and q^. The d!^. operators and are connected by momentum conservation, id!^^ + 
idt^^ = q^, so one of them can be eliminated. Hence we expect that reparametrization 

invariance will give a large number of connections on the homogeneous basis, so many in 
fact that it is not even convenient to write down the homogeneous basis. It is much quicker 
to just write only the RPI basis and expand it to determine a basis of allowed homogeneous 
operators. 

The RPI basis for three jets at LO is made of two quarks fields and a gluon field (we 
do not consider here the case with pure gluon jets), rii and n2 will be the directions of the 
quark and antiquark jets, and will be the direction of the gluon jet. As for the two jet 
case, because of current conservation, the only objects that can carry the vector index and 
are RPI invariant are gl^ and 7^. The RPI basis is 

Jl = ^n^,Co^lvillqaigQnl,^3^n2,Co2 , ^2 = ^ n^^Co^ll^ Ivq J gQnl,Co;'^ n2,U2 , 

J3 = ^ ni,C,^l<TlulligGnl,Cjz^ n2,Cc2 J4 = ^ ni,C,^lllulaigGnl,CJz^ 5 



J9 = ^,.^ni,<ii^«^^n3,<i3^^12<7^n2,<i2 > JlO = ^"1,^17^!^ 7/^^^^^^3,i3^^12<7*n2,<i2 , 

Jll = ^n„C^,llq.tgQZ,^-M^^-^r^2,^2 ■ (120) 

For the first four operators we chose the Dirac structures {luilj^qa, I'^il^qa, lalulj^, 
ijilula} in order to simplify the transformation of the basis under charge conjugation. Since 
{7J, ^} = the sum of the first two structures gives —gj^^, iq^, and using the antisymmetry 
of Q'Z the sum of the last two gives AgJ^j^j^^, so structures with a gj^^ are redundant. Other 
three 7 operators are also redundant. We have used the equations of motion and Bianchi 
identity in Eqs. (15111541) to eliminate and momentum conservation to eliminate id^^ = 
q^ — id^_^ — idf^^. For the operators Jg-n we have a derivative contracted with igG^s, 
and we can use the gluon equation of motion, idn^^Gns — i^o- — "^•^mo- — '^dn2a)Qns = • • •) 
where the ellipsis denotes higher twist terms, to eliminate {id^^ + id^^) leave only 
id 12 = (^^ ni ~^'9n2)- Note that we cannot use the trick used in DIS for Q^, to eliminate 
Jg-ii, because here id^_^^ and have different collinear directions. Operators with two or 
more derivatives are redundant for the construction of the LO basis of RPI currents with 
one-vector index /i, and hence do not need to be considered. 

We can match the three jet RPI basis of currents with the basis of homogeneous SCET 



37 



currents by writing 



r^^ ^„ ■J A r^, r.^^ ■J A 



ni,n2,n3 



ni, ".2,713 



;i2i) 



On the RHS the integration variable was changed using d^i = n-quji and any additional n-q 
factors were absorbed into the Wilson coefficients Cf {(jJi). We can determine the currents 
[j7/(t^j)]3_jet of the homogeneous basis, whose form is as in Eq. (l86l) . by just expanding the 
currents fll20p using Eqs. (IHTj) and (!83l) . This yields the homogeneous operator basis 



Ji 

J2 



It Xn2,uj2 ' 

\n2,^2 1 

J-i = ^3 Xn^,LoA9^nz)^3i3lTXn2,u^2 > 
J4: = ^3 Xni,uiilTi''i{9^nz)'-'3Xn2,L02 i 

— Xni,LL)i1^iT{.9^n-i)^^-iXn2,u)2 i 

— ^2 Xrn,u)i_f^2T^9^nz)^zXn2,ijJ2 i 

Jl = ^1 Xni,a;i<r^1^3(fl'^n3)a;3Xn2,a.2 , 
Js = ^2 Xm,cuini^T^i'-ii9^nz)uJsXn2,iU2 > 



(122) 



n3Ti9l3, 



'n3±/'^3 



(n2j.u;2+^ii.u;i)Xn2, 



,<.^2 5 



Ju Xni,u)i 7r(5'^n3±)<.^3("'2j/^^2 + 'T-lj/^l)Xn2,tJ2 ) 



n2,uJ2 1 



where n'^rp 



rij* — q^ {riyq) / q"^ , n^rp = — q^{n2-q)/q^ and ^ = (^3 ■ 5)^3/2 + (ns ■ 9)^3/2. To 
simplify the results we did not bother to write out the terms with qfi{9B!^^j_) in Eq. (11221) . 
which are terms that vanish in a frame where g_L„3 = 0. In some cases we have absorbed 
RPI factors in the Wilson coefficients Ci{(jJi) when carrying out the expansion. 

The tree level matching from QCD to SCET for three jets comes from matching two 
Feynman diagrams in QCD onto the operator basis in Eq. (I122p . and is done at the hard 
scale n = Q. This gives 



Ci — Cq 



-2 



-C, 



-2 



n2-n3 UJ2UJ2, 



a 



3,4 



C7- 



11 



0. 



(123) 



The results for these Wilson coefficients are invariant under type-III RPI as expected. 

The above matching computation can be compared with the tree level SCET compu- 
tations for parton showers in Ref. 57|], where three final state jets are considered. To 
compare the calculations we take the two stages of matching of Ref. 57|] both at fi = Q, 
and we split the operators in Eqs. (27,28) of Ref. [57| into two parts, O3 = O^a + C'sb and 



(2) 



•-^30 ^ ^ 



and O'f}^ and a frame qn^_ 

+ C^J^ = O^a 
results in Eq. ([1231) . 



(2) 
3b • 



The matching computation of Ref. 57| used a frame q. 



m2- 



for O 



3a 



(2) 

for O^h and . With these frame choices, we con&m that 



Ci + C,J, = Osa + of} and C2J2 + C5 J5 



Osb + Og^ , providing a cross-check on the 
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G. Two Jets from Gluon Fusion: gg qq operators 

Next we consider the example of the production of two quark jets from gluon fusion, which 
is relevant for the LHC In this application we will see that RPI substantially constrains the 
number and structure of operators. This basis of operators have not yet been constructed. 



The factorization theorem for pp 2 jets has been discussed in Ref. |58|, and were also 



considered recently in Ref. |59| using SCET. SCET has also been used to resum electroweak 



Sudakov logarithms by solving RGE equations for four quark coUinear operators in Refs. 60, 



61l . |62| | , and to consider Higgs production from pp collisions [63 . 

We consider the incoming gluons to be coUinear in different directions, which is appro- 
priate for the high energy collision of energetic protons at the LHC, and we assume that the 
final state jets have a large perpendicular momentum relative to the beam axis. Hence the 
final jets are described by two additional collinear directions, making four in total. Unlike 
our previous examples, here there is not an external vector, the hard interaction takes 
place entirely between strongly interacting particles. Hence this is an example of the case 
ii) discussed above Eq. (I23l) . 

Similarly to the three jets case, it is convenient to directly write the RPI basis without 
first writing the homogeneous basis, because the presence of four collinear directions imply 
that there are a large number of homogeneous operators, many of which are restricted by 
RPI. Due to the absence of an external hard vector g'^ in this process, in the definition of the 
currents we make use the RPI delta function factors of Eq. (13T|) . Akm- The general formula 
for matching the RPI operators onto homogeneous operators is 

ni,n2,n3,n4 i i,j km ni,n2,n3,n4 i i 

+ ..., (124) 

where we use the same manipulations needed to get Eq. (127|) . Note that here we have divided 
the RPI operators into the 5-functions in Akm which depend on uJkm, and the remainder of 
the operator that does not. The starting point for building a basis for is the object 
^ni^^3 ^nf ^"2- We assume a LL + RR chirality for the quarks which is suitable when 
strong interactions produce massless quarks, and hence include either 7''* or 7'^7°"7'^. Since 
the overall operator is a scalar, all the vector indices on the field strengths and on the Dirac 
structure must be contracted with g^^^s or idl^^^s. We can use the equations of motion and 
Bianchi identity in Eqs. fl51ll53ll54p to eliminate terms with in any operator, and terms 
with dn^^Gns dru^iQn^^- In addition, momentum conservation implies idif^_^ + idi^^ + idf^^ + 
i(9^^ = 0, and we will use this to eliminate all operators with an idm- This leaves twenty 
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operators for the RPI basis 



Qi 




Q2 


= ^mlf.g.a^dn.ptgGii'^igg^^^n^ , 


Qs 




Q4 


= ^mlf3guaigGii''Jggn!idn2f,'^n2 , 


Qs 




Qe 


= ^n^lf,lylJdn:,pigQZ^9Gf^^n2 , 


Qt 




Qs 


= ^n,lulalpigQZigQn^,^dn2t.^n2 , 


Q9 








Qio 


= ^mlf.idnsaigGnsigOn^idn^uidn^/s'^n^ ■ 




(125) 



The other ten operators Qii„2o have the same structure as Eq. (11251) but with a trace 
over color for the gluon operators, for example Q^^^ = ^ni'~^i3gua^'i^[igGnl^dniiJiigQni]^n2- 
Note that Qi_io have Qnz to the left of ^„4, so one might think that there are ten more 
operators with the ^'s in the other order. However, in Eq. (11241) we sum over 113^4 and 
integrate over duj^duj^, and hence include operators obtained from the interchange ^ 
n^, uj^ <-> UJ4. Recall that the directions are only determined by the matrix elements. 
So if we consider a matrix element with gluons in the n and n' direction then there is a 
contribution from n^ = n, n4 = n', and from n^ = n\ 71,4 = «. Other possible operators might 

be ^ni'yaidn3pigQii^idni^,igGn4^9n^'^^n2, ^m1t^idn3aigQii^idn4uigGniidn2(3'^n2 and similarly 
with the trace. We can use the Bianchi identity (JMl) to rule them out. For example, in the 
first operator we have implicitly already used the Bianchi identity for the idn^^^igQ!^^ term 
because we did not write operators with ic^migQn^- But we can apply the Bianchi identity 
to idn^pigQt^'^ , that is not connected with 7's. In this way we can write this operator in 
terms of Q^, Q4 and operators with three gluon fields. Note that we do not need to consider 
operators with idn ■ idn' since all these contracted derivatives are contained in the Akm^s. 

A natural frame for analyzing gg — >■ qq is the CM frame with the choices ni=n2, ^^2=^1, 
77-3 = 77-4, 714 = 123. We expand the currents (I125P with an eye towards using them in this frame. 
Actually, only the condition n3 = n4, 114 = 713 is necessary to find the following operators 

(126) 

U2 = Xni,uJi is {gi3'^s±)co:Agi3ni^) ul4Xn2,u)2 ' 

03 = LO2 Xni,uJi{g^n3i.)uii{gn2-B^^)t^^Xn2,ui2 i 

04 = ^2 Xn^,uj^{gn2-B^^)uj.j,{gt^nil.)uiiXn2,U2 ! 
^5 = ^iXnx,0Jiil'A{g^n3±)LO3{g^nA±)uiAXn2,u)2 ? 

= ^?,Xni,oJiil''i{g^n3±)u)3{g^nA±)oJ4Xn2,u)2 5 
Oj = UJ2UJZUJ4 Xn^^u^jkii {gl$mi)ujz{gn2-B^^)^^Xn2,W2 , 
Os = l^2^:i^^iXn^,u>jkii{gn2-B^^)^.Xgl$nil.)uAXn2,uJ2 ) 
O9 = {uJ2fuJ4Xrn,ujrii ig7l2- B^^)^^{g7l2- B^J^^Xn2,u^2 , 
OlO = {uJ2f(^:iXrn,ujri^ ign2- B^^)^^{g7l2- B^J^^Xn2,u^2 ■ 

Oii_2o have the same structure of (I126p but with a trace over color of the two gluon operators. 
Oi is given by the expansion of Qj for i = 1,2,5,6, by the expansion of a suitable linear 
combination of Qj and Qj_i for z = 3,6, and of Qj and Qj_3 for i = 4,8. O9/10 are given by 
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the expansion of a suitable linear combination of Q9/10, Q1/2 and Q4/3. In some cases we 
have absorbed reparametrization invariant prefactors that appear in the expansion into the 
Wilson coefficients Ci{uJi). By using momentum conservation it is possible to reduce these 
ten operators to just four independent operators at leading order in SCET.^ 

It is straightforward to carry out the matching from QCD onto the SCET operators in 
Eq. (11261) . At tree level there are three Feynman diagrams. The amplitude squared is also 



known analytically at one-loop [6J] , and a full matching computation at this order involves 
regulating infrared singularities in the same way for the loops in QCD and SCET before 
subtracting. The only point to be careful about is the sum over the nj's in Eq. (I124p . since 
definite values for these n^'s should be determined by the states. For example, if we consider 
the tree level gg qq matrix element of Oi with perpendicular polarization for the gluons 
then 



+ ig^C,{uj[, uj',)u;', Kl^e^^.^^,] [u^,^ri',T^T^u^,^] . (127) 

The two terms come from the cases 71,3^4 = n'^ 4 and n3^4 = '^^4 3 respectively. Therefore to 
determine the C^'s it suffices to compute terms contributing to the color structure T'^T^ in 
QCD, which at tree level gives 

Cl = 7 ^ , C2 = 7 ^ , C3 = , C; 



(n3-n4)to'3tU4 ' (ra3-n4)co'3U;4 ' (^2-^4) to'2t^4 ' (^2 " ^4) U;2to'4 ' 

C4 = C6_2o = . (128) 

Note that the results for the C^s are invariant under type-Ill RPI transformations as ex- 
pected, and that in the frame used for our computation 77,3 ■ 77,4 = 2. We have confirmed 
that a consistent result is obtained by considering the T^T^ terms. Eq. (I127p expresses the 
interesting fact that with distinct coUinear directions for all final state particles, only the 
color ordered QCD amplitudes are needed for the matching which determines the SCET 
Wilson coefficients. 



V. CONCLUSION 

In SCET the momenta of collinear particles are decomposed with light-like vectors n'^ 
and n'^, where n is close to the direction of motion. The vectors n'^ and n'^ are required to 
define collinear operators that have a definite order in the power counting. However, there 
is a freedom in defining n and n, which leads to reparametrization constraints. The decom- 
position of operators in the theory must satisfy these constraints in order to be consistent. 
This reparametrization invariance gives nontrivial relations among the Wilson coefficients 
of collinear operators occurring at different orders in the power counting, and for situations 



^ We thank W. Waalewijn for his expHcit derivation of this point. 
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with multiple collinear directions gives constraints on the form of operators making up a 
complete basis. 

In this paper we have constructed objects that are invariant under both reparametrization 
transformations and collinear gauge transformations, a superfields for fermions and a 
superfield Q^'^ for gluons. Here the subscript Ui denotes an equivalence class of light-like vec- 
tors under RPI. The superfields are invariant under collinear gauge transformations through 
a reparametrization invariant Wilson line Wn, that is the generalization of the usual Wn^ ■ We 
constructed RPI operators out of these superfields by introducing reparametrization invari- 
ant 5-functions. The 5-functions act on the RPI operators to pick out large momenta, and 
are convoluted with hard Wilson coefficients that must be computed by matching computa- 
tions. The power of the RPI operators is that they encode information about the minimal 
basis of Wilson coefficients. However, they do not have a definite power counting order. By 
expanding them in A one obtains a minimal basis of operators with a good power counting, 
where all constraints on the Wilson coefficients are made explicit. The final basis of opera- 
tors with a good power counting involves a two-component field Xn^ for quarks, a field 
for the two physical gluon polarizations, derivatives and delta functions 6{u — 

that pick out the large momenta of these collinear fields. That is it. Other field components 
such as rii-Bni, and other derivatives such as irii-dm, are eliminated from the purely collinear 
operator basis using the equations of motion. 

This procedure was applied to several processes. We studied spin-averaged DIS for quarks 
at twist-4, as a means of testing our setup in a framework where the power suppressed 
basis of operators is well understood. We then constructed a minimum basis of pure glue 
operators for DIS at twist-4. These applications involve a single collinear direction. For 
processes with multiple collinear directions we considered operator bases for jet production. 
Useful constraints from RPI were not found for the first power suppressed operators in 
e~^e~ 2 jets. On the other hand, already at leading order in the power counting, RPI 
provided important constraints on the complete basis of operators for e~^e~ — 3 jets with 
three distinct collinear directions. RPI was also very useful in constructing a complete basis 
of operators for gluon fusion producing two quark initiated jets, where there are four collinear 
directions. In this case the process of interest is pp 2 jets, which will be studied at the 
LHC. We expect the complete bases of operators constructed here will be a useful ingredient 
in the study of factorization theorems for this process. The steps we used to construct 
complete basis will also be useful when considering factorization for processes with more 
jets in the final state. In general we found that RPI becomes more powerful for processes 
involving more jets, essentially because the number of vectors and proliferates faster 
than the number of objects that must be considered to build the RPI basis. 

An interesting observation discussed in section IIVGI is that when matching from QCD 
onto SCET operators describing multiple collinear directions rij, the Wilson coefficient is 
determined by the color ordered QCD amplitude. Since results for multi-leg QCD amplitudes 
are often expressed in a color ordered form, this should simplify the matching of QCD 
amplitudes onto SCET. 

This work was supported in part by the Director, Office of Science, Office of Nuclear 
Physics of the U.S. Department of Energy under the Contract DE-FG02-94ER40818. IS 
was also supported in part by the DOE OJI program and by the Sloan Foundation. 



42 



APPENDIX A: INVARIANCE TO THE CHOICE OF HARD- VECTOR q^' 



From the construction in section IIII A[ a natural question arises about the special role 
of in Eq. flSTl) . What happens if there is more than one possible choice for q^ in a given 
process? Say we have a q^ and a q'^ with Wilson coefficients that can depend on g^, g'^, and 
g ■ g', where g^ ~ g^ ~ A. It turns out that in this case any linear combination of q^ and q'^ 
in Eq. (ETj) is equally good, and is equivalent to any other choice. Hence, one choice suffices. 
To prove this we consider the expansion of the reparametrization invariant variable 



W (2^y_4^^^^ . (Al) 
g2 \ ^ q^ q n ■ q 

where we take the plus sign if the expansion is done with n ■ q' /n ■ q > n ■ q' /n ■ q and the 
minus sign otherwise. Now use this variable to define 

q' ■^^-i{q■^^)^Qfl^, (A2) 
where the operator qInv R-P-'- ^"^^ expansion starts at order A^. Thus 

m 



duo C{uj) 5{uj — q' ■ id) = J duj C{lo) 5[uj — ^ q ■ id — QInv) 
= I diu' Ci^u') 6iuj' -q-id- qISv/0 

= j duj' C{uj') 5{uj' -q-id) + B{lu') Q[3v 5'{uj' -q-id) + ... , (A3) 

where in the second line we changed the dummy variable to u' = ^u, In the last line both 
terms are RPI, and the ellipsis denotes higher order terms which are also RPI order by order 
in A. Eq. flA3P demonstrates that we can swap the parameter g' — > g in the 5-function, since 
the change is compensated by a change of notation in the leading order Wilson coefficient 
C ^ C . Given that we imagine starting with a complete basis of RPI operators built with 
5{uj — q'-id) or with 6{uj' — q-id), the higher terms in the series in Eq. flA3p . like B, also 
simply change a Wilson coefficient in our basis. Thus, the choice of g or g' in the 5-function 
just corresponds to a different choice of the basis for the invariant operators, and one choice 
suffices. 



[1] J. C. Collins, D. E. Soper, and G. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988), 

hep-ph/0409313. 
[2] G. Sterman (1995), |heFph/9606312. 



[3] C. W. Bauer, S. Fleming, and M. E. Luke, Phys. Rev. D63, 014006 (2001), hep-ph/0005275 



[4] C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart, Phys. Rev. D63, 114020 (2001), 
hep-ph/0011336, 



[5] C. W. Bauer and I. W. Stewart, Phys. Lett. B516, 134 (2001), hep-ph/0107001 



[6] C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D65, 054022 (2002), hep-ph/0109045| 



43 



[7] C. W. Bauer, S. Fleming, D. Pirjol, I. Z. Rothstein, and I. W. Stewart, Phys. Rev. D66, 
014017 (2002), hep-ph/0202088: 

J. Chay and C. Kim, Phys. Rev. D65, 114016 (2002), 'hep-ph/0201 1971. 

M. Beneke, A. P. Chapovsky, M. Diehl, and T. Feldmann, Nucl. Phys. B643, 431 (2002), 



hep-ph/0206152 



R. J. Hih and M. Neubert, Nucl. Phys. B657, 229 (2003), hep-ph/0211018 



C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D68, 034021 (2003), hep-ph/ 0303156 

M. E. Luke and A. V. Manohar, Phys. Lett. B286, 348 (1992), hep-ph/9205228, 

A. V. Manohar, T. Mehen, D. Pirjol, and L W. Stewart, Phys. Lett. B539, 59 (2002), 

hep-ph/0204229l 



D. Pirjol and I. W. Stewart, Phys. Rev. D67, 094005 (2003), hep-ph/0211251 



C. M. Arnesen, J. Kundu, and I. W. Stewart, Phys. Rev. D72, 114002 (2005), hep-ph/0508214[ 
A. Hardmeier, E. Lunghi, D. Pirjol, and D. Wyler, Nucl. Phys. B682, 150 (2004), 
^hep-ph/0307171. 



J. Chay, C. Kim, and A. K. Leibovich, Phys. Rev. D72, 014010 (2005), hep-ph/0505030 



R. J. Hih, T. Becher, S. J. Lee, and M. Neubert, JHEP 07, 081 (2004), [hei>ph/0404217 
H. D. Politzer, Nucl. Phys. B172, 349 (1980). 
R. L. Jaffe and M. Soldate, Phys. Lett. B105, 467 (1981). 
R. L. Jaffe and M. Soldate, Phys. Rev. D26, 49 (1982). 
R. K. Ellis, W. Furmanski, and R. Petronzio, Nucl. Phys. B207, 1 (1982). 
R. K. Ellis, W. Furmanski, and R. Petronzio, Nucl. Phys. B212, 29 (1983). 
R. L. Jaffe (1996), hep-ph/9602236, 

F. J. Tackmann, Phys. Rev. D72, 034036 (2005), hep-ph/0503095. 



A. V. Manohar and I. W. Stewart, Phys. Rev. D76, 074002 (2007), |hep^ph70605001 , 
C. W. Bauer, C. Lee, A. V. Manohar, and M. B. Wise, Phys. Rev. D70, 034014 (2004), 
hep-ph/0309'278l 



J. Chay, c7Kim, Y. G. Kim, and J.-P. Lee, Phys. Rev. D71, 056001 (2005), [hep^ ph/0412110| 
C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D67, 071502 (2003), 'hep^ph/0211069j 
M. Beneke and T. Feldmann, Phys. Lett. B553, 267 (2003), hep-ph/0211358, 
M. Neubert, Phys. Rept. 245, 259 (1994), hep-ph/9306320, 



S. Fleming, A. K. Leibovich, and T. Mehen, Phys. Rev. D68, 094011 (2003), |hei>ph/0306139 
L Z. Rothstein, Phys. Rev. D70, 054024 (2004), [hep^ih70301240. 



A. K. Leibovich, Z. Ligeti, and M. B. Wise, Phys. Lett. B564, 231 (2003), 'hep-ph/0303099' 
M. Finkemeier, H. Georgi, and M. Mclrvin, Phys. Rev. D55, 6933 (1997), hep-ph/9701243. 
C. L. Y. Lee (1997), hep-ph/9709238. 

B. Geyer and M. Lazar, Phys. Rev. D63, 094003 (2001), hep-ph/0009309. 

A. Retey and J. A. M. Vermaseren, Nucl. Phys. B604, 281 (2001), hep-ph/0007294. 

S. J. Brodsky, D. S. Hwang, and I. Schmidt, Phys. Lett. B530, 99 (2002), hep-ph/0201296. 

J. C. Collins, Phys. Lett. B536, 43 (2002), hep-ph/0204004, 

X.-d. Ji, J.-p. Ma, and F. Yuan, Phys. Rev. D71, 034005 (2005), [hep^ph/0404183[ 
A. V. Belitsky and A. V. Radyushkin, Phys. Rept. 418, 1 (2005), hep-ph/0504030. 
J. B. Kogut and D. E. Soper, Phys. Rev. Dl, 2901 (1970). 
S. Bashinsky and R. L. Jaffe, Nucl. Phys. B536, 303 (1998), hep-ph/9804397. 
J. Bartels, C. Bontus, and H. Spiesberger (1999), hep-ph/9908411. 

C. W. Bauer, A. V. Manohar, and M. B. Wise, Phys. Rev. Lett. 91, 122001 (2003), 



hep-ph/0212255 



44 



C. Lee and G. Sterman, Phys. Rev. D75, 014022 (2007), [hep^h/0611061 

C. Lee and G. Sterman (2006), hep-ph/0603066f 

M. Trott, Phys. Rev. D75, 054011 (2007), hep-ph/06083001 

S. Fleming, A. H. Hoang, S. Mantry, and L W. Stewart, Phys. Rev. D77, 114003 (2008), 
10711.20791 

M. D. Schwartz, Phys. Rev. D77, 014026 (2008), 0709.2709. 

S. Fleming, A. H. Hoang, S. Mantry, and L W. Stewart, Phys. Rev. D77, 074010 (2008), 
hep-ph/0703207. 

T. Becher and M. D. Schwartz (2008), 0803.0342. 

A. Jain, I. Scimemi, and I. W. Stewart, Phys. Rev. D77, 094008 (2008), 0801.0743. 
A. H. Hoang and S. Kluth (2008), 10806.38521 

C. W. Bauer and M. D. Schwartz, Phys. Rev. Lett. 97, 142001 (2006), hep-ph/0604065. 
C. W. Bauer and M. D. Schwartz, Phys. Rev. D76, 074004 (2007), hep-ph/0607296. 
N. Kidonakis, G. Oderda, and G. Sterman, Nucl. Phys. B525, 299 (1998), hep-ph/9801268. 
C. W. Bauer, A. Hornig, and F. J. Tackmann (2008), 0808.2191. 

J.-y. Chiu, F. Golf, R. Kelley, and A. V. Manohar, Phys. Rev. D77, 053004 (2008), 0712.0396. 
J.-y. Chiu, F. Golf, R. Kelley, and A. V. Manohar, Phys. Rev. Lett. 100, 021802 (2008), 
0709.2377. 

J.-y. Chiu, R. Kelley, and A. V. Manohar, Phys. Rev. D78, 073006 (2008), 0806.1240. 
V. Ahrens, T. Becher, M. Neubert, and L. L. Yang (2008), 0808.3008. 
R. K. Ellis and J. C. Sexton, Nucl. Phys. B269, 445 (1986). 



45 



